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SPIRAL STELLAR NEBULAE AND COSMIC GAS JETS* 


BY 


C. E. R. BRUCE! 


SUMMARY 


The paper describes the successful culmination of a search for photographic 
evidence of the existence of spiral stellar nebulae, initiated by the writer's electrical 
discharge theory of astrophysical atmospheric phenomena, as applied to some types 
of bright line stars. This theory for the first time offers a satisfactory explanation 
for the general form and the gas movements observed in planetary nebulae. It 
elucidates the general nature of cosmic gas jets, concerning which quite erroneous 
assumptions are made in relation to those in the sun’s atmosphere which cause mag- 
netic storms. Applied to galactic nebulae, the theory explains the strange discon- 
tinuity recently observed in the gas movements in the barred spiral nebula NGC 1097. 


1, SPIRAL STELLAR NEBULAE 


The forces which are operative in the formation of the extra-galactic 
spiral nebulae have remained a mystery since the first discovery (1)? of 
their spiral structure, be it noted, by Lord Rosse, President of the Royal 
Society in 1845, using his 6-ft. diameter “leviathan of Parsonstown.”’ 
To these the writer has recently added for like consideration spiral 
nebulae on a stellar scale (2a), the existence of which had remained 
apparently unnoticed and without comment since they were photo- 
graphed at Lick Observatory by H. D. Curtis (3) in 1918. 

The existence of these stellar spiral nebulae was recently ‘‘discov- 
ered”’ as the result of a search of the literature for just such photographic 
evidence of their existence. 

The first step which originated this search was the suggestion that 
the well-known group of P Cygni stars, which have been the special 


* Based on several recent Electric Research Association Reports, including Ref. Z/T 117, 
Z/T 125 and Z/T 132. 

1 Electrical Research Association Laboratories, Leatherhead, Surrey, England. 

* The boldface numbers in parentheses refer to the references appended to this paper. 

(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
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study, for example, of C. S. Beals, Dominion Astronomer of Canada, 
are stars essentially similar to R Monocerotis, better known as Hubble’s 
Variable Nebula, shown in Fig. 1. This photograph has the additional 
interest that it was the first official photograph taken by Hubble him- 
self using the 200-in. Hale telescope at Mount Palomar Observatory on 
January 26, 1949. However, in the P Cygni stars the associated come- 
tary nebula, it is suggested, is directed toward us. 

On the electrical discharge or “‘thunderstorm’’ theory of these 
astrophysical atmospheric phenomena (2b, c) these are stars in which 
the atmospheric gas has been largely collected along the discharge 
channel of an old nova-like outburst by magnetic pinch effect, just as 


Fic. 1. Hubble’s variable nebula (R Monocerotis). 


the pressure is increased in the laboratory arc discharge by the same 
effect. It is therefore of note that P Cygni is indeed an old nova. It 
is this magnetic pinch effect which, for example, mainly renders possible 
electric arc welding. The pressure gradient so generated from the 
current constriction at the cathode spot along the axis of the arc, 
causes a jet of hot gas to flow and entrain the metal of the welding elec- 
trode and carry it over to the anode, as will be seen in Fig. 2, a photo- 
graph taken from a high speed film of the welding process taken by 
Needham in these laboratories. 

This axial pressure rise in an electrical discharge is proportional to 
the product of the current and the current density. Both these quan- 
tities decrease as the discharge is propagated outwards in these stellar 
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atmospheric electrical discharges, so that outward jets of hot gas emit 
the bright lines in the stars’ spectra which are therefore displaced to- 
wards the violet in these stars in which the nebulae are directed towards 
us, as it is suggested it is in P Cygni. 

There exist, however, many stars in which the bright emission 
lines are double, with components displaced towards both the red and 
the violet. These, it seemed to the writer, could be stars in which there 
had been two major discharges in the original outburst, just as is always 
Ghseeved in the galactic spiral nebulae (2d). Since photecraphs of 
stars like Hubble’s Variable Nebula exist, and will account for these 
stars having a spectrum of single bright emission components, it seemed 


Fic. 2. High speed photograph of the arc welding process. 


quite likely that some stars might sometime have been photographed 
with two such old discharge channels, if the proposed explanation of this 
type of stellar spectrum is correct. Enquiry from a few astronomical 
friends brought from Dr. Merrill at Mount Wilson Observatory a refer- 
ence to the paper by Curtis already referred to (3), which contains a 
profusion of examples of just such pairs of old discharge channels as the 
theory had suggested might exist. Some of these are shown in Fig. 3, 
and it was quite exciting to find that in some of these the two discharge 
channels had assumed a spiral form, for example NGC 4361, while in 
others like NGC 2452 and NGC 6778, they showed a barred-spiral form, 
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exactly analogous to what is observed when the corresponding outbursts 
occur on a galactic scale. The view of spiral galactic nebulae has been 


elaborated in a recent E.R.A. Report (2d). In contrast to the galactic 
discharges, the two discharge channels remain straight in the majority 


N.G.C. 2392 


“NGC. 2371-2 


2452 


N.G.C. 6853 


Fic. 3. Photographs of planetary nebulae (Curtis (3)). 


of stellar outbursts. In others they remain straight over a proportion 
of their length before turning off into a spiral form, when the result is a 
barred-spiral nebula. 

This duality of discharge channel is such a marked characteristic of 
these planetary nebulae that Curtis himself had offered an explanation 
in terms of matter falling in towards the star from the two poles. It 
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will be seen that, in contrast to the explanation now offered, this would 
lead to the red component of the lines being the brighter in those stars 
possessing a spectrum of pairs of bright components, since it would be 
emitted by the arm on the star’s nearer side. Dr. Merrill in two 
papers (4a, b) has referred to twenty stars possessing spectra with 


Figures along axés are the N 

velocities in km/sec of the 

Rand V components of the — | 

spectrum lines 101 
‘ 


"135 
23 + 126 —— 
26 


Full lines and©--indicate discharge jets coming towards us. 
Dotted lines and+ indiccte jets directed away from us. 


Fic. 4. Diagrammatic representation of the gas jets in NGC 2392. 


pairs of bright emission components in which he could compare the 
relative brightnesses of the two components. In thirteen of these the 
violet component was the brighter, in six the two components were 
approximately equal, and in only one was the red component of the 
lines the brighter, in accord with the prediction of the new theory, and 
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quite the reverse of what is to be expected were Curtis’s suggestion 


correct. 
2. GAS MOVEMENTS IN THE PLANETARY NEBULAE 


The theory also affords an explanation for the gas movements ob- 
served in these nebulae, a study of which was made at Lick Observatory 


3 


Fic. 5. Moore’s drawings of the N; line in NGC 2392, with the slit along 
(a) the minor axis and (b) the major axis. 


by Campbell and Moore (5) and published at the same time as Curtis's 
photographs. The movements in many of the nebulae defied all 
efforts at their explanation on the basis of any of the many hypotheses 
tried by Campbell and Moore, but are, like the photographic and spec- 
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troscopic observations, also in accord with the new theory. In the 
diagram in Fig. 4 the figures along the axes are the gas velocities given 
by Campbell and Moore for NGC 2392, the type nebula of the group 
which particularly resisted all attempts at their explanation. These 
will be seen to be at once explicable as deriving from the two gas jets 
generated by the two diametrically opposed discharge channels. The 
high velocities of the two jets are seen superposed in the central regions, 
so that the spectrum lines are there split as in Fig. 5. It is only in the 
SE and NW quadrants that the two jets are separable. In the diagram 
in Fig. 6 the figures in the circles are average values derived from Camp- 


NW 


N 


= 


SE 


Figures in circles are gas velocities in km/sec. 


Fic. 6. Diagrammatic representation of the gas jets in NGC 6543. 


bell and Moore’s measurements for NGC 6543 obtained with the slit 
of the spectrograph in ten positions across the image of the nebula. 
These show that in this star the discharges gave rise to spiral jets of 
gas, as the photograph of the nebula also indicates. 


3. COSMIC GAS JETS 


So far as the writer is aware, the magnetic pinch effect in electrical 
discharges is the only physical process which can gather atmospheric 
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Fic. 7. Photograph of the lines Ha and [N I] ] 46583 in the spectrum of the barred spiral 
nebula NGC 1097. Reproduced with permission of the publisher, from the Astrophysical 
Journal, Vol. 132, p. 30, copyright (1960) by the University of Chicago Press (8). 


matter into the discrete arms found in both the stellar and galactic 
spiral nebulae, and at the same time give rise to the discrete jets of gas 
observed in both these types of nebulae, and in the solar atmosphere, 
where they account for magnetic storms. The existence of these solar 
jets has been postulated to account for terrestrial magnetic storms since 
1896, but their origin has hitherto remained a mystery, and quite 
erroneous suppositions have been made, as we shall see, regarding their 
physical characteristics, suppositions which the writer believes are 
physically impossible. For it is quite erroneous to assume, as has often 
been done in the past, particularly in these theoretical studies of mag- 
netic storms, for example, that the temperature and velocity of these 
gas jets can be treated as entirely independent characteristics. 

On the contrary, the velocity of propagation of the jet can only 
derive from the velocity of sound in the gas forming the jet, and it 
seemed to the writer that this relationship can be used as the basis of a 
cosmic ‘‘gas-velocity thermometer’’ (2e). This ‘thermometer’ has 
been used to ‘‘measure’’ temperatures over the range 5000° K. to 
4,000,000,000° K., and its ‘‘readings’’ will be seen in Table I (2f) to 
agree well with other methods of estimating the temperature. 


TaBLE I.—Gas Temperatures and Velocities in Cosmical Electrical Discharges. 


Temperature, ° K. Velocity, km./sec. 
Source Theory Observation Theory Observation 
Long-period —_ 5000° to 10,000° 8.5 to 12 9 (Joy) 
variable 11 (Merrill) 
stars 
Combination 110 (Merrill) 
spectra star 8 X 10°° (Fe VII, Fe X (?)) —_— 
(AX Persei) 
Extra-galactic 108° ~ 109° _ 1750 to 5400 1800 to 4250 (Seyfert) 


nebulae 
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Indeed another of the notable recent predictions of the theory was 
based on this ‘‘gas-velocity thermometer,’’ as was recently mentioned 
in a note in ‘Reader Opinion’’ (2g). It derived from the observed 
velocity of the solar jets which enter the upper atmosphere, and which 


Fic. 8. Photograph of the barred spiral nebula NGC 1097, taken at McDonald Observa- 
tory. Reproduced with permission of the publisher, from the Astrophysical Journal, Vol. 132, 
p. 30, copyright (1960) by the University of Chicago Press (8). 
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cause magnetic storms and aurorae. These velocities, according to 
Meinel (6), reach 3500 km per second. 

It was pointed out in a contribution to the Institution of Electrical 
Engineers Convention on Thermonuclear Processes in April last year 
that these velocities indicate that the temperature in these solar dis- 
charges must reach over 100,000,000° K. between the sun and the 
Earth. This prediction had not long to await observational confirma- 
tion, for satellite observations made by scientists of the U. S. Navy last 
year showed that solar flares are accompanied by the emission of 80,000 
V X-rays, indicating, as they point out, the existence of temperatures 


of 100,000,000° K. (7). 


#, THE GAS JETS IN BARRED-SPIRAL NEBULA NGC 1097 


Support for the theory of spiral and barred-spiral nebulae put 
forward in the above sections, as applied to galaxies, is contained in a 
recent paper (8) by Burbidge and Burbidge on the ‘‘Motions in Barred 
Spiral Galaxies.”” The authors show a photograph, reproduced in 
Fig. 7, of the emission lines Ha and [N II] \6583, taken with the slit 
of the spectrograph more or less along the bar of the barred-spiral 
galaxy NGC 1097 (Fig. 8). 

These emission lines are seen each to comprise two displaced por- 
tions separated by the nucleus. The authors’ Table 1 of the gas veloci- 
ties observed along each of these lines shows that these velocities are 
constant along the bar out to distances ranging from 2’’.9 to 15’’.8 on 
each side of the nucleus, but that the velocities in the line of sight on 
the two sides of the nucleus differ by about 415 km per sec. As the 
authors say “there is essentially a velocity discontinuity in the inner 
part.” 

This will be seen to be in complete accord with the view of stellar 
and galactic barred spiral nebulae discussed in the present paper. The 
two electrical discharges give rise to diametrically opposed jets of gas, 
flowing outwards along the bar from the center. It is these two jets of 
gas which give rise to the two displaced portions of the spectrum lines. 
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ELLIPTIC FUNCTIONS IN NETWORK SYNTHESIS 


BY 


S. D. BEDROSIAN'! 


ABSTRACT 


Admittedly the Jacobian elliptic functions are not too familiar to the engineer. 
However, there are occasions when alternate techniques are not as direct or as con- 
venient to solve the approximation problem of network synthesis. The role of 
elliptic functions in providing a common thread for fundamental concepts involved 
in Cauer image parameter filters, filters on an insertion loss basis and wide band 90° 
phase-difference networks is shown. A “unified” design chart for these networks is 
also included. Additional tabular and graphical data are presented emphasizing 
extreme parameter values not found in the literature, that is, modular angle @ > 89°, 
to facilitate getting numerical results. An attempt has been made to provide a prac- 
tical basis for engineering insight into the synthesis problem. It is also hoped that 
this paper may serve as an inducement to some of the engineers interested in network 
synthesis to acquire broader acquaintance with these topics via the references cited. 


1, INTRODUCTION 


Functional relations between the parameters of Jacobian elliptic 
functions are well known (1—4).?. Since the pioneering work of Legendre 
many of these have been tabulated because of their importance in 
applied mathematics (4-12). The advent of digital computers has 
made significant extensions possible (13, 14). In spite of this, their 
applications to the approximation problem of network synthesis have 
been handicapped by the coarseness and inadequacy of tabulations 
hitherto considered to be detailed and complete. The difficulty stems 
from the doubly periodic nature of the elliptic functions which generally 
necessitates double entry tables (15—16). 

In such network applications the values of the elliptic functions of 
particular interest and difficulty are those for the modulus k approaching 
unity. (See Appendix I for a brief summary.) A review of the rela- 
tionships between parameters, the methods of tabulation and the tabu- 
lations themselves reveal that adequate means are at our disposal to 
surmount many of the handicaps by a judicious choice in presenting 
the results. Table I and the curves in Figs. 1, 2 and 3 show the rela- 
tionships germane to network synthesis. In some cases it will be 
feasible to get the necessary data by graphical interpolation of such 
curves drawn to much finer scale. Often it will be sufficient to use 
existing tables augmented by just a few additional calculations with 
the formulas provided. 

Although a fairly lengthy bibliography has been included, it is by 
no means exhaustive. Rather it is hoped that it will serve as a point 


1 The Institute for Cooperative Research, University of Pennsylvania, Philadelphia, Pa. 
* The boldface numbers in parentheses refer to the references appended to this paper. 
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of departure for the several topics touched upon in this concise treat- 
ment of a broad subject. 


2. APPROXIMATION PROCEDURE 


The approximation problem of network synthesis is concerned with 
the derivation of a transfer function based on performance specifications. 
An approximation is usually required since the specifications tend to be 
idealized and often are not quite realizable as a finite physical network. 
Considerable ingenuity and mathematical knowledge have been brought 
to bear in this area. 

The application of Jacobian elliptic functions to the approximation 
problem originated with Cauer in 1931 (18). He addressed himself 
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Fic. 3. Elliptic function parameters for extreme values of modular angle. 


to the problem of achieving an equiripple, or Tchebycheff, approxima- 
tion to constant attenuation in the stop band of image parameter filters. 
Subsequently (1939), they were used to obtain equiripple character- 
istics in both the pass and stop bands of filters synthesized by the inser- 
tion loss theory (19). By then, however, modular transformations 
were introduced on the elliptic functions to simplify the resulting ex- 
pressions. These elliptic function approximation methods were also 
applied (1950) in the design of wide band phase-difference networks 
(20, 21). Winkler gives a comprehensive survey of the status of this 
and other approximation techniques as of 1954 (17). Note also the 
methods given by Papoulis and Deutsch (31, 32). 
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The key point to observe, for our purposes, is that the elliptic 
function filters and the phase-difference networks have three distinct 
frequency bands of interest. The former has equiripple attenuation in 
the pass and stop bands separated by a monotonic transition region. 
The latter has the inverse situation, with a prescribed equiripple phase 
characteristic over a given bandwidth, on either side of which there is 
a monotonic phase region. This is shown in Fig. 5. 

We should note here that the filter discussion is being confined to 
the low pass case. This involves no loss of generality inasmuch as the 
usual transformations can be applied to get the high pass, band pass 
and band stop cases if desired. 


3. DESIGN PARAMETERS 


The pertinent elliptic function parameters have already been dis- 
cussed fairly extensively in connection with the specific synthesis pro- 
cedures. Useful supplementary material has also been provided by 
Glowatski (22), Grossman (23), Henderson (24) and this author (25). 
As general references, there are the several more recent texts on network 
synthesis and the special ‘‘Modern Filter Design Techniques Issue”’ 
of the JRE Transactions on Circuit Theory (26-30). The latter serve 
to bring the reader up to date through 1958. 


A. Image Parameter Filters 


It must be emphasized that the similarity between the insertion loss 
filters and the Cauer image parameter filters is limited to the theory 
involved in the choice of design parameters. The discrepancy in actual 
insertion loss characteristics is due to interaction effects in the pass 
band and reflection effects in the rejection band. Nevertheless the 
Cauer method did represent a significant improvement over the con- 
ventional Zobel image parameter filter design technique in that it 
permits independent consideration of the image impedance and the 
propagation functions. Network realization, although usually in 
ladder form, is based on the symmetrical lattice structure. 

The image impedance for the symmetrical lattice is 


Z1(w) = + VZ.(w)Z,(w). (1) 


The transfer function, P;(w), can be expressed in terms of a parametric 
ratio function given by 


(2) 
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== Ar(w) 
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or alternately as 


H(w) = tan 


(3) 


The independent parametric relations of Eqs. 1 and 2 determine a 
purely reactive symmetric lattice. The positive roots are normally 
used (26). 

It is well known that for a low pass filter Z;(w) is controlled by the 
cutoff frequency w, and the stop band critical frequencies, while H(w) 
depends on the cutoff and pass band critical frequencies. It is con- 
venient to use a normalized frequency parameter x = w/w. The 
properties of the transfer loss A;(x) are intimately related to those of 
H(x). In fact the stop band frequencies at which H(x) is unity de- 
termine the poles of A;(x), that is, the frequencies of infinite attenua- 
tion. The maximum possible number of poles for A1(x) is one greater 
than the number of pass band critical frequencies. Thus specification 
of these poles and cutoff frequency determines A;(x) as well as H(x). 
(See Fig. 4.) 
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The information required to start a Cauer parameter filter design 
is tolerance e, with respect to the minimum value of attenuation amin 
versus the corresponding coverage k. That is the frequency at which 
Qmin is first reached. (See Fig. 4.) Then we have 

2 


(4) 


‘) 
Qmin = In 
€ 
2 


For close tolerance, that is, «small, these can be simplified to « = Te 


or its equivalent 


and amin = In- respectively. From this it is necessary to determine 
€ 


the number and location of the control frequencies in the pass. band. 

Location of the control frequencies for Tchebycheff approximation 
(or geometric mean variation with respect to unity) to a flat transfer 
loss is facilitated by a change in the independent variable in the H (x) 
function of the form x = k sn u, where u = F(¢, k); also, x = ksin ¢. 
Then the range of the frequency variable 0 < x < k is equivalent to 
0<@<7/2 or 0<snu <1 or O < F(g,k) < K, as shown in 
Appendix I. 

Then we have the pass band control frequencies 


x ; = sn 
A~An+1-—-i n + 1 
and poles of A;(x) in the stop band 


1 
Xp(2(n—i)+3) = (1/k) sn [ (22 1)K/2(n + 1)] 


1<i<n+1 (6) 


Once the number and location of the control frequencies of A;(x) are 
known, the flat loss characteristic of the filter may be computed (26). 

A similar treatment applies to Z;(x) where the design specifies a 
tolerance on the deviation of the impedance function Z(x) = Z;(x)/R 
from unity over a specified range of frequency. It also can be consid- 
ered a special case of the function F:, shown in Fig. 4. The number and 
location of the control frequencies in the stop band are required. In 
this case the stop band control frequencies are given by 
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and the unit values of Z(x) on the normalized frequency scale x are 
given by 


2(n —1) +3 


= ksn 


B. Elliptic Function Filters 


On the other hand, for the elliptic function filters, based on insertion 
loss theory, we start off by letting the required transfer function be a 
general rational function, given by 


1 
|Zi2|? = 1+ (9) 


Then the desired filter is specified by any three of the following four 
parameters: the minimum loss in the stop band, the maximum attenua- 
tion in the pass band, the width of the transition region, and the com- 
plexity of the function n. This is based on the selection of the ap- 
propriate Tchebycheff rational function R,(w) which is of the form 


R.(o).= (1 — (1 — (1 
or (10) 
R w (ws? w?) athe (Wen? _ w?) 
= (1 — wow?) (1 — wa2w*)- (1 — 


where the first is even and the second is odd. Each function has zeros 
on the w axis in the range 0 < w < 1 and poles in the complementary 
range 1 < w < ©, the zeros and poles, respectively, being located at 
reciprocal points on this axis. Because of the reciprocal relationship 
between the zeros and poles, it is clear that R,(1/w) = 1/R,(w). Note 
that here 7 is an integer which determines the complexity of the function. 

The number of ripples in the pass band and the number in the stop 
band are equal and are determined by the number of poles n. The 
transition interval is a function of the selectivity of the filter as shown 
in Fig. 5. The transition interval becomes narrower as functions of 
higher order m are used. Of course for finite networks this interval 
cannot be reduced to zero. It is convenient to define a reference fre- 
quency, w) = V¥w,w., which is somewhat analogous to the cutoff fre- 
quency used in image parameter filter design. Again, as a matter of 
convenience, we work with a normalized frequency variable 2 = w/wo 
= w/Vw,w.. Then, as indicated in Fig. 5, the pass band terminates 
at @ = Vk, and the attentuation band begins at @ = 1/Vk, where 
k = w,/w, is the selectivity parameter, by definition. 
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The corresponding insertion power ratio required for a symmetrical 
filter turns out to be of the form 


ee = 1 + — 1)[R(Q) (11) 
where 
R(2) = 8S 


Loss 


Insertion 


73 


Phase Difference 


Fic. 5. Typical elliptic function filter and phase-difference network performance 
curves showing three distinct frequency bands. 
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The following additional requirements must also be met: 
(a) |R(Q)| = — 1) (e» — 1) = 1/ky 
(b) So = 1/Vk, 


(c) The selectivity parameter k and the discrimination parameter k, 
are related by the condition (2M + 1)K’/K = Ki'/K,. An alterna- 
tive expression for this condition is given in terms of modular constants 


as qi = g@™+, 
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(d) The frequencies at which the loss is zero and the reciprocal of the 
frequencies at which the loss is infinite are given by 


s=0,1,2,---,M 


2M +1 n=2M +41. 


In the usual case, the filter design problem consists of the specification 
of the number of ladder filter sections M required to attain a prescribed 
minimum band loss aq, a prescribed transition band vk <2 < 1/Vvk, 
and a tolerable pass band ripple a,. The well known expression re- 
lating the primary design parameters is 


aq = 10[log (e%» — 1) — (2M + 1) log gq] — 12.04 db. (13) 


For additional special design formulas see (19, 23). 


C. 90° Phase-Difference Networks 


The third type of network involves derivation of two transfer im- 
pedances, suitable for all-pass networks, and having approximately 2/2 
radians difference between their phase angles over a prescribed band of 
frequencies. These are usually called 90° phase-difference, or phase- 
splitting, networks. It can be shown that a suitable rational function 
where poles and zeros are the negatives of one another can be written as 


tan (6/2) = (14) 


7=1 + 
where 


= 1/2 — — B:). 


This function can be written in terms of parametric equations to give a 
Tchebycheff approximation to a constant 90° phase difference. Normal- 
izing the prescribed band of frequencies so that the geometric mean of 


the ends of the band Vw4w,z equals unity, the parametric equations are: 


tan (6/2) = Vk, cd [(2nK,/K)u; ki] (15) 
and 
w = 1/Vvk’ dn [u; (16) 


If we express the maximum value of phase difference @ by ®,,, we can 
rewrite Eq. 15 as 


tan (#,,/2) = (15a) 
Since k; «< 1 normally, we can apply the approximation 
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The corresponding modular constants are related by 


qi = gq". 


Then Eq. 15 assumes the alternate form 


tan (,,/2) = 2q". 


Practical applications are confined to small error angles so that 


®,,/2 = tan (#,,/2). 


These approximations result in an equation relating phase tolerance, 
bandwidth and network complexity. 


®,, = 49” (17) 


or alternately 


= 4[exp (— 1K’/K) (17a) 


where 


®,, = phase-tolerance angle in radians 
n = network complexity, that is, the number of first order all-pass 
sections in network N plus the number in network P. (See 
Fig. 6.) 
g = modular constant or “‘nome”’ related to the bandwidth of the 


network by sec 6 = 1/k’ = we/wa. 


Applying some transformations of the elliptic functions to the Eq. 16 
we get the necessary equation for determination of the value of the 
poles and zeros required for all-pass phase-difference network synthesis. 


cn(u;; k) 
wn/wa) sn(u;; k) 


= (wg/wa)'cs(u;; k) 


and 


j 


The theoretical sideband rejection attainable by use of these networks 
can be expressed in terms of a given error angle ®,, as 


(19) 


Rejection (db) = 20 log [cot (#,,/2) ]. 
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4. PARAMETER DETERMINATION 


In the previous section we summarized the already known funda- 
mental design parameters of the three types of networks being consid- 
ered. In this section are presented means of relating these parameters 
involving elliptic functions. This is done in a way which facilitates 
determination of numerical values while at the same time emphasizing 
the inherent similarities in concept of these types of networks. 


EC Case for n=6 


Network 
La 


Network 


Fic. 6. Phase-difference network configurations. 


Specialized representations have been devised by various authors to 
show the performance capabilities of these networks for various ranges 
of parameter values including network complexity. Because of the 
common thread of elliptic functions utilized, however, it becomes 
possible to use a single representation for all three types at the expense 
of three sets of scales. This is shown in Fig. 1. If desired the ranges 
given therein can readily be extended. The alternate representations 
are also amenable to similar modification. 

Similar simplification and consolidation of information are possible 
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in connection with location of the critical frequencies which are neces- 
sary for synthesis of the networks. A family of curves typified by 
Fig. 2 practically permits determination of these critical frequencies by 
inspection. An added feature here is that the curves go considerably 
beyond the usual limitation of modular angle @ = 89°. A tabulation of 
“amplitude” @ to twelve decimal places is available in the Spenceley 
tables for values of r and @ in integral degrees (10). 

The basic nature of the curves shown in Fig. 2 can better be ap- 
preciated when it is realized that, at least in principle, it is possible to 
get from them the Cauer parameters (given in the Glowatski tables (22)) 
and also the normalized solutions for pole-zero pairs for all-pass phase- 
difference networks (given by Bedrosian (25)). 

To get the parameters given in the Glowatski table for the plot, first 
refer to Table III to get the specified values of r for the desired n from 1 
through 12. Next read off the values of sin ¢ from Fig. 2 for these 
specified r and the desired modular angle 6. Now we have only to 
multiply these readings of sin¢@ by vk for the appropriate 6, where 
k = sin @, to get 


= Vsin 6-sn(mK/n; k) m=1,2,---,n—1. (20) 


Note again that Fig. 2 permits extension of the tables. 

To obtain the normalized solutions for phase-difference networks 
from the plot, first refer to Table II to get the specified values of r for 
the desired » from 1 through 12. Next read off the values of cot ¢ 
from Fig. 2 for these specified r and the desired modular angle 6. (The 
bandwidth ratio 1/k’ is related to 6.) Now we have only to multiply 
these readings of cot ¢ by vk’ for the appropriate 6, where k’ = cos 8, 
to get 


xj =es{ — j = 0,1,2, — 1. (21) 


A few words are in order regarding the source of Tables II, III, 
and IV. Observe that the elliptic functions being utilized are of the 
form sn(u;k), cs(u;k), etc. Furthermore, the variable u is normally 
given in terms of the complete elliptic integral K which, depending 
on the elliptic function in question, is a half period or a quarter period 
of the function. To facilitate the use of the Smithsonian Elliptic 
Function Table note that the Spenceleys let 


u = (r/90)K (22) 


and tabulate it as r (in degrees). In other words the independent vari- 
able u = F(¢, k) is divided into 90° of a quadrant K, that is, into 90 
equal steps. In addition to tabulating the usual elliptic functions for 
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TABLE IV.—Tabulation of Variable u in Terms of Angle r in Degrees for Image 
Parameter Filter A pplication. 


(Where N + 1 = n, the number of unit values of H(x) or Z(x) and N is the 


number of control frequencies.) 
r 


45.000 

30.000 60.000 Use with Fig. 2 and eqs. 5 and 7 

22.500 45.000 67.500 (Note that they are reciprocal) 

18.000 36.000 54.000 72.000 

15.000 30.000 45.000 60.000 75.000 

12.857 25.714 38.571 51.428 64.286 77.143 

11.250 22.500 33.750 45.000 56.250 67.500 78.750 


N w/2(N+1) 24/2(N+1) 34/2(N+1) 44/2(N+1) 5/2(N+1) 62/2(N+1) 72/2(N+1) 


N 

1 Use with Fig. 2 and Eqs. 6 and 8 

2 45.000 (Note that they are basically reciprocal except 
3 33.750 56.250 for ordering of values obtained) 

4 S 27.000 45.000 63.000 

5 : 22.500 37.500 52.500 67.500 
6 
7 
N 


19.286 32.143 45.000 57.857 70.715 
16.875 28.125 39.375 50.625 61.875 73.125 


a/4(N+1) 54/4(N+1) 77/4(N+1) 92/4(N+1) 114/4(N+1) 137/4(N+1) 


r = 0 tor = 90 for each degree of the modular angle @ (from 6 = 1 to 
6 = 89), they list the corresponding value of ¢. This ¢ is the amplitude 
of u. Thus the values of r corresponding to the values of u in the 
pertinent equations are listed in Tables II, III, and IV. 

Finally we consider Fig. 3 which, in conjunction with Table I, 
permits comparisons of various combinations of the parameters K, 8, k, 
k’, g, etc. as found in the numerous tabulations cited in the Introduction. 
Here again the range of values presented exceeds that found in the 
literature. These results should aid the researcher requiring the use of 
these elliptic parameters at the high values of 6 —- 90° and k > .999. 
For example, the range of the modular constant is 0 < gq < 1 and of the 
complete elliptic integrals is 7/2 < K < «. The values given here 
go up to g > 0.8 and K > 24. With the choice of scales in Fig. 3, 
it would be a simple matter to extend the plots as required and still 
retain the ability to use effective graphical interpolation. 


5. CONCLUSIONS 


The illustrations given in the preceding sections should be sufficient 
to demonstrate the important role of elliptic functions in the initial 
approximation step of network synthesis. The field of application of 
the tabular and graphical data furnished is by no means limited to the 
three types of networks discussed. The necessary additional back- 
ground and detailed discussion of synthesis procedures are beyond the 
scope of this paper. The insight provided should facilitate understand- 
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ing and use of the material available in the references. Particular 
attention is invited to the increasing number of tables of normalized 
element values useful in meeting practical design problems (30). 


APPENDIX I 


Elliptic function formulas: 


dy 


i Panty Elliptic integral of the first kind (1-1) 


u = F(g,k) = 
dx 
u = F(x, k) 


where 1 > k > 0 is the modulus, ¢ is the amplitude of the elliptic integral (that is, am(u, k)), 
and x = sing. When ¢ = 7/2 in Eq. I-1: 


(1-3) 


(1-2) 


Similarly if x = 1 in Eq. 1-2: 


dx 
Fi,k) = K = (1-4) 


There are three basic elliptic functions: 

sn(u,k) = sing =x Odd function (1-5) 
(1-6) 
(1-7) 


cn(u, k) = = vi — 


functions 


dn(u,k) = Ad = V1 — Rx? 


Note that these functions are trigonometric functions of ¢, the amplitude of u. When con- 
sidered as functions of u they are doubly periodic and are called elliptic functions. The inde- 
pendent variable u is not an angle. It may be real or imaginary as shown in Table V. 


TABLE V.—Periods of the Basic Functions. 


Function Real Period Imaginary Period 


sn(u, k) 2jK’ 
cn(u, k) 2(K + jK’) 
dn(u, k) 4jK’ 


The analogy to trigonometric functions is further shown by the relations: 
sn?u + cn®u = 1 


dn?u + k*sn?u = 1. 


| 
| 
Also 
i 
| 
| 5 
: 
(1-8) 
| 
(1-9) 
i 
k sau 
k 
cn u 7 : 
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Figure 7 shows the interrelations. Recall the limits on the modulus are 
We now define the complementary modulus k’ 


Also (see Fig. 8) 
k = sin@ 


k' = cos0 


where 6 = modular angle. (N.B. x/2 — 0 = @ —a in Fig. 8.) Then 


dy 
Fi(¢, k’) = 
J, v1 — (k’)? sin? y 


When ¢ = x/2 
F(x/2,k') = K’. 


Next we define g, the modular constant or “‘nome”’ 
q = exp(— rK’/K). (1-15) 


This basic number g may also be obtained from the modular angle 6. Graphical representa- 
tions of these functions for real u are determinable in a simple way from the plots of the elliptic 
integral as a function of the amplitude. Using Fig. 9, we can get the results shown in Table VI 
for @ = 60°: 

TABLE VI. 


sn(u, k) cn(u, k) 


:1>k>0. 
(1-10) 
(1-11) 
(1-12) 
| 
ee 0 0 0 1 1 oe 
K r/2 1 0 k' =} 
3K —1 0 =} 
4K de 0 1 1 = 
= 
“2 
or}. 
a=! 
° 
gis 
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APPENDIX II 


To get numerical values for the complete elliptic integrals of the first kind K = F(x/2, k) 
and K’ = F(x/2, k’), we offer an alternate form of the series expansion based on the Wallis 
Theorem (2, 11): 


=5[1+2(F )+9(%) +50(F y+ +2208(F) +---] for 1. 


For @ > 89.95°, it is simpler to use the expression given by Kaplan (9) to get results with 
ten-place accuracy: 


Kin (1+Z)ing - 


Numerical values of the modular constant g can be obtained by use of the Weierstrass ap- 
proximation 
q = d + 2d® + 15d® + 150d3 + 1707d"7 20,910d + 268,616d%* + --- 


also 


where 


Vcos 1 — vk’ 
@ 1+ vk’ 


The Spenceley tables (10) indicate that this approximation is not good for sixteen-place ac- 
curacy for @ > 70°. For larger modular angles we can use: 


Inging’ = 7 
or alternately 


log log ¢’ log log = 1,8615228349, 


where q’ is the modular constant of the complementary modulus k’. 
if 


REFERENCES 


(1) E. T. WittaKER ANp G. N. Watson, “Modern Analysis,” 4th ed., New York, Cambridge 
University Press, 1940, pp. 429-535. 

{ (2) H. Hancock, “Elliptic Integrals,” republished by Dover Publications Inc., 1958. 

(3) A. Eacue, ‘The Elliptic Functions As They Should Be,” Cambridge, England, Galloway 
and Porter, 1958. 

(4) A. M. LeGEnpRE, ‘‘Traite des fonctions elliptique,’’ Tome 2 Paris 1826, pp. 222-243. 
Tables reprinted in 1931, by F. Emde Stuttgart, and in 1934 by the Biometrika Office, 
London. 

(5) K. Hayasut, ‘“‘Taflen der Besselechen, Theta, Kugel und anderer Funktionen,” J. Springer, 
Berlin, 1930. Table of K, pp. 72-81. 

(6) L. M. Mitne-Tuompson, “Ten Figure Table of the Complete Elliptic Integrals,’ Proc. 
London Math. Soc., Vol. 33, pp. 160-164 (1930), reprinted by Dover Publications as 
“Jacobian Elliptic Function Tables,” 1950. 

(7) J. R. Arrey, “Toroidal Functions and the Complete Elliptic Integrals,” Phil. Mag. 
Ser. 7, Vol. 19, pp. 177-188 (1935). 

(8) C. HeumAN, “Tables of Complete Elliptic Integrals,” J. Math. Phys., Vol. 20, pp. 127-206 
(1941). 

(9) E. L. Kapvan, ‘Auxiliary Table of Complete Elliptic Integrals,” J. Math. Phys., Vol. 25, 

pp. 26-36 (1946). 


| 
| 
i 
if 
‘ 
i 
: 
: 
} 


30 S. D. (J. 


(10) G. SPENCELEY AND R. SPENCELEY, “Smithsonian Elliptic Function Tables,’ Wash. D. C., 
Smithsonian Institution, Vol. 109, 1947. 

(11) E. JAHNKE AND F. Empe, “Tables of Functions with Formulae and Curves,” 4th Ed., 
reprinted by Dover Publications, Inc., 1945. 

(12) M. ScHULER AND H. GEBELEIN, “Acht und neunstellige Tabellen zu den elliptischen 
Funktionen,” Berlin, J. Springer, 1955. 

(13) A. R. DiDonato anp A. V. HerRsHEY, ‘‘New Formulae for Computing Incomplete 
Elliptic Integrals of the First and Second Kind,’’ NPG Report No. 1618, NAVORD 
Report No. 5906; 30 Jan. 1959. 

(14) R. G. SELFRIDGE AND J. E. MAxFieELp, “A Table of the Incomplete Elliptic Integral of 
the Third Kind,”” New York, Dover Publications Inc., 1959. 

(15) H. J. Orcnarpn, “Computation of Elliptic Functions of Rational Fractions of a Quarter- 
period,” IRE Trans. on Circuit Theory, Vol. CT-5, pp. 352-355 (1958). 

(16) L. Gotp, ‘“‘New Approach to Elliptic Functions,’’ JouR. FRANKLIN INstT., Vol. 267, 
pp. 503-520 (1959). 

(17) S. WINKLER, ‘‘The Approximation Problem of Network Synthesis,”” JRE Trans. on Cir- 
cuit Theory, Vol. CT-1, pp. 5-20 (1954). 

(18) W. Cauer, ‘“‘Theorie der linearen Wechselstromschaltungen,”’ Berlin, Akademie-Verlag, 
Ger., 2nd ed., 1954. Also in English, ‘Synthesis of Linear Communication Networks,” 
New York, McGraw-Hill Book Co., 1958. 

(19) S. DarLincton, “Synthesis of Reactance 4-Poles,”’ J. Math. Phys., Vol. 18, pp. 257-353 
(1939). 

(20) S. DARLINGTON, “‘Realization of a Constant Phase Difference,’’ Bell Sys. Tech. J., Vol. 29, 
pp. 94-104 (1950). 

(21) H. J. Orcnarp, ‘Synthesis of Wideband Two-phase Networks,” Wireless Engineer, 
Vol. 27, pp. 72-81 (1950). 

(22) E. GLtowatsk1, “Sechstellige Tafel der Cauer-Parameter,"’ Munich, Verlag der Bayer- 
ischen Akademie der Wissenchaften, 1955. . 

(23) A. J. GRossMAN, ‘Synthesis of Tchebycheff Parameter Symmetrical Filters,’ Proc. TRE, 
Vol. 45, pp. 454-473 (1957). 

(24) K. W. HENDERSON, ‘“Nomographs for Designing Elliptic Function Filters,’’ Proc. TRE, 
Vol. 46, pp. 1860-1864 (1958). 

(25) S. D. BepRostAn, ‘‘Normalized Design of 90° Phase-Difference Networks,"’ JRE Trans. 
on Circuit Theory, Vol. CT-7, pp. 128-136 (1960). 

(26) M. B. REEp, ‘Electric Network Synthesis,’’ Englewood Cliffs, Prentice-Hall, Inc., 1955. 

(27) E. A. GuILLEMIN, “Synthesis of Passive Networks,"’ New York, John Wiley & Sons, Inc., 
1957. 

(28) J. E. Storer, ‘‘Passive Network Synthesis,’"” New York, McGraw-Hill Book Co., 1957. 

(29) D. F. Tutte, Jr., ‘Network Synthesis,” Vol. I, New York, John Wiley & Sons, Inc., 
1958, Chap. 14. 

(30) Modern Filter Design Techniques Issue, JRE Trans. on Circuit Theory, Vol. CT-5, (1958). 

(31) A. Papou.is, “On the Approximation Problem in Filter Design,”” JRE National Con- 
vention Record, Vol. 5, Part 2, pp. 175-185 (1957). 

(32) S. Deutscn, ‘‘A Different Approach to the Approximation Problem,”’ Proc. IRE, Vol. 48, 
pp. 1175-1176 (1960). 


x 
4 
4 
: 
r 
: 
: 


A METHOD OF DETERMINING PHASE CENTERS AND ITS 
APPLICATION TO ELECTROMAGNETIC HORNS* 


BY 
¥. ¥. HU! 


ABSTRACT 


A method of locating the phase center of any radiating system from the expres- 
sion of its radiating field is formulated. This method is then applied to electro- 
magnetic horns of different dimensions and flare angles. It is believed that the 
results and discussions presented in this paper will be useful in the design and position- 
ing of the feeding horn such that the paraboloidal reflector will produce a desirable 


radiation pattern. 
1, INTRODUCTION 


One of the important applications of electromagnetic horns is the 
use as primary feeds of paraboloidal reflectors. In the usual practice, 
the relative locations of the horn with respect to the reflector are ad- 
justed by a cut-and-try process until the measured directivity of the 


radiation pattern is “‘maximum.” It is also known in practice that a 
“maximum” directive radiation pattern can be obtained if the field 
distribution on the aperture is equiphase. By the geometry of the 
paraboloidal reflector, equiphase aperture illumination can be achieved 
by locating a point source as the primary feed at the focal point. It is 
apparent that an electromagnetic horn is not a point source, in fact, 
there is no point source in nature. However, any radiating system is 
equivalent to a point source within a certain angular region if the equi- 
phase surfaces of its radiation field are concentric spheres in that region. 
The center of these concentric equiphase spheres is defined as the phase 
center of this radiating system. In this paper, a method of locating 
the phase center of any radiating system from the expression of its 
radiation field is formulated. This method is then applied to electro- 
magnetic horns of different dimensions and flare angles. It is believed 
that the results presented in this paper will be useful in the design and 
positioning of the feeding horn such that the paraboloidal reflector will 
produce a desirable radiation pattern. It is further believed that the 
method of locating phase center developed here will also be useful in 
determining the phase center of any other radiator when such informa- 


tion is desired. 


* This work was supported in part by Rome Air Development Center, Rome, N. Y., 


in 1954 under Contract No. AF30(602)-926. 
1 Electrical Engineering Department, Syracuse University, Syracuse, N. Y. 
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Il. A METHOD OF LOCATING PHASE CENTERS 


A method of locating the phase center of any radiating system was 
developed by the author during a study of the phase center of electro- 
magnetic horns (1).2- For the convenience of presentation, a conical 
horn is used as an illustration. Referring to Fig. 1, let the coordinate 
origin be the vertex of the conical horn. If there exists a phase center 
of the horn, it must be located on the horn axis due tosymmetry. The 
distance from the phase center to the origin is denoted by d which is to 
be determined. Let the spherical coordinates of any field point P be 


Pr. 6.4) 


Fic. 1. Phase center and coordinate system. 


(r, 8, @) with 6 = 0° as the horn axis. By the cosine law, the distance R 
between the field point P and the phase center can be written as 


(1) 


The expression is independent upon the azimuth angle @. In the far 
field, d <r, Eq. 1 can be approximated by 


R = [r + — 2rd cos 


R=r-—dcos@. (2) 


The approximate relation (2) can also be obtained by considering 
distances r and R as parallel lines when P is in the far field. 
The field with such a phase center can be expressed as 


(3) 


where | M| is the amplitude of the field. The phase angle » depends also 
upon the choice of time reference. Substituting (2) and k = 22/) into 
Eq. 3, we have 


E, = 


E, = | (4) 


2 The boldface numbers in parentheses refer to the references appended to this paper. 
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On the other hand, from the theory of electromagnetic radiation, the 
far field produced by a finite size radiating system can always be ex- 
pressed in the following form 


E, = |M(r,0, (5) 


If the radiating system has a ¢-symmetry, that is, the radiation field 
E, is independent of ¢, then Eq. 5 can be written simply as 


By comparing the two phase factors in Eqs. 4 and 6, it can be con- 
cluded that: if the field expressed by Eq. 6 is spherical, then h(@), 
generally, contains two terms, one is (dcos 6)/\ and the other is a 
constant 6. The term (d cos 6)/\ indicates the spherical nature of the 
wavefront. The constant term 6 depends upon the choice of the time 
reference. By suitable choice of the time reference, 5 can be made 
equal to zero. Inthe case é = 0, the ratio h(6)/cos 6 will give the value 
of d/\ which determines the location of the phase center. A more 
general way of obtaining d/) is to plot h(@) versus cos 6; it should be a 
straight line passing through the origin for 6 = 0 and its slope is the 
value d/X. When h(@) contains an additional non-zero value of 4, 
the plot method can still be used. In such a case, the straight line will 
not pass through the origin ; however, its slope still gives the required d/h. 

Note that the plot of 1(@) versus cos 0 is a straight line over the entire 
range of 6 (0° < 6 < 180°) only if the radiator is an ideal point source. 
For a practical radiating system the plot is a straight line only within a 
certain angular region. Within that angular region, the radiation field 
is spherical and the location of the phase center, d/X, is the slope of the 
straight line part. If no part of the plot is a straight line, then no part 
of the wave front is spherical, and the concept of phase center is mean- 
ingless in such a case. 

It is apparent that this method of determining the phase center is 
very general and can readily be applied to any radiating system with 
its radiating field expression known. In this paper, the application of 
this method to sectoral horns (1) only is presented. Besides sectoral 
horns, this method was later also used to determine the phase centers 
of helical beam antennas (2,3). 

For a radiating system without a ¢-symmetry, a single phase center 
is not defined. However, for any ¢ = constant plane, the phase center 
can be defined and determined in the same way. The wavefronts in 
such a plane are simply circular arcs. For sectoral horns, the phase 
centers in the two principal planes are of practical interest and will 
be treated in this paper. 
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Ill, RADIATION FIELDS OF SECTORAL HORNS 


Rectangular horns have been studied in considerable detail in the 
literature (4-7). The expressions of the radiation fields of sectoral 
horns have been worked out by Barrow and Chu (4). Their results 
are derived by using the aperture-diffraction method. In the deriva- 
tion, the aperture surface is taken to coincide with the cylindrical wave- 
front of the dominant mode at the mouth of the horn. The aperture 
field is assumed to be that of the incident wave with the effect of the 
reflected wave neglected, and the radiation field is then computed by 
applying the aperture diffraction formula. The field expressions for 
the £ and H sectoral horns in the two principal planes, @¢ = 0° and 
@ = 90°, are listed below in the forms which are convenient for the 
present applications. 


Fic. 2. E-plane sectoral horn. 


1. E-Plane Sectoral Horn (Fig. 2) 
(a) in the H-plane, ¢ = 0° (abbreviated as EH case) 


cos (= sin é 
E, = 4+(2)(*) f (Br) cos 6 
T r a? —ag 


1 — 45, sin’ 6 
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(6) in the E-plane, @ = 90° (abbreviated as EE case) 


2 
+ V1 (>) (Bre) cos (6 cos (6-0 19’, (8) 


Fic. 3. H-plane sectoral horn. 


2. H-plane Sectoral Horn (Fig. 3) 
(a) in the H-plane, ¢ = 0° (abbreviated as HH case) 


E, = Es (277) (cos p6’) cos (0 — 6’) 


+ jH,?” (25 (cos cos (0-8) 79’. (9) 


(b) in the E-plane, ¢ = 90° (abbreviated as HE case) 


sin sin a) 
E, = Aj Ain ikr (2) (272? ) (cos p60’) cos 6’ 


2x r 


+ (2 cos (cos po) | 


where p = 1/2a) and H, is the Hankel’s function of the second kind 
and pt“ order. The meaning of the remaining symbols appears in Fig. 2 


or Fig. 3. 


| 

| 

| 

E = 

7 

tase | 

| 

of 

i 

a 

: 

| 

| 

620 

| 

r 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

y 

| 

| 

| 


*(aueld @ SOD snsaaa 


*(auel @ S09 snsaaa 


b 


| 


T 


WP 
£0990 -00000/ _ | 9810/ 

| -LEL50 2 800 - 


oo, 


| | 


ines 
| 


2099. 0| 626280 | | 
- + 060260 £8860 | 
$/16/80\ \ GE | | 


‘642090 | 


“£0999 0 E/PBE/ | 


1€9060 | €2068/ Tes 


696 66 0 | 10494 } - 

259660 

£65 9 | 


98266 0 
66660 
a0007 | 
| 


| 


161502100 | 


20 O02 


- + + + + — — 
| 
bi 
| | 18] SE 
| | 
| 3 8 | | 
8 
y 
2 e | | 
2 4 | | 
| | | 
| | | | | | 
| 
= ~ ~ ~ Ss 
—— 


*(auerd 7777) @ S09 snsaaa *(aueld FH) @ S09 snsz90 


007 066: 096" oor 26° 


UP £69 =P / 


0- 0 PP SEP O -20/03/ 


802 - 9802 206 


/ 


| £20660 - 00000) = 7 £0997 0 -00000/ 


/ 


/ avoid FH 
6L209 | 
669700 
PEFR O £0999 0 | PPSEP 
180860 \ 1€9060\ 
420660 696€6'0\ 
299460 665960\ 
9SL660\ 199860 1 
6€6660| 619660\ O12¢¢% 7 
000001 2010s 0 000001 2g10S7 
@ 502 907 


n 
= 
< 
a 
= 


i. 


ob! 


Jan., 1961. 
» 
ee 
' 
. 
[ 
: 
s Ss 
(0/4 
= 


38 Y. Y. Hu (J. F. 1. 


All the integrals in Eqs. 7, 8, 9 and 10 are complex functions of a 
real variable 6, and can be written in the following polar form. 


Integral = | A (6)|e?"*. 


Each integral combined with e-**" and the other factors can always be 
written in the same form as given by Eq. 6 which is the desired form 
for locating the phase center. 


IV. RESULTS AND DISCUSSIONS 


Phase centers have been computed for sectoral horns of different 
geometrical dimensions. Four sample plots of h(@) versus cos 6, one 
for each case, are presented in Figs. 4 to 7 to show how each phase 
center is determined. Four tables of computed 4(6@) and cos @ are listed 
on the figures for reference. In each plot, a straight line can be drawn 
from 6 = 0° up to a certain angle. The straight line relation means 
that the wavefront in that range of @ is circular and the slope of the 
straight line gives the location of the phase center d/d. It is noticed 
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Fic. 8. Phase center location as a Fic. 9. Phase center location as a 
function of the flare angle (E plane sectoral function of the flare angle (H plane sectoral 
horn). horn). 


that in the two flared plane cases, EE (Fig. 4) and HH (Fig. 7), the 
straight line covers only a small range of 8, about — 8° < @ < 8°; 
while in the two non-flared plane cases, EH (Fig. 5) and HE (Fig. 6) 
the ranges of @ are much wider. This is connected to the fact that the 
main beam in the flared plane is narrower than that in the non-flared 
plane. 

Compvted phase centers for E and H sectoral horns are presented in 
Figs. 8 and 9 to show how the location of the phase center varies with the 
flare angle a. These phase centers are computed for a/d = 0.7, 
r2/\ = 5 and a from 5° to 20°. In both figures, the following proper- 
ties are observed. 


1. At small flare angles, the phase centers in the E-plane and H-plane 
almost coincide. They deviate from each other as the flare angle 
increases. 
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2. In the non-flared planes (EH and HE), the phase center remains 
almost fixed and is located close to the aperture of the horn even when 
the flare angle varies in the other plane. 

3. In the flared planes (EE and HH), the phase center moves inward 
as the flare angle increases. 


From the above properties, some suggestions can be made about 
how to use the phase center information of sectional horns in producing 
approximately uniform phase illumination on the aperture of a para- 
boloidal reflector. 


1. For horns with phase centers in the two principal planes located 
close to each other, satisfactory operation can be obtained by locating 
the focus of the reflector at the midpoint between the two phase centers. 

2. For horns with phase centers in the two principal planes located 
not too close to each other, satisfactory operation can be obtained by 
locating the focus at a point between the two phase centers such that 
the maximum phase deviations in the two planes are approximately the 
same. This point can be determined by the use of the geometric 
optics method. 

3. For sectoral horns, as a rough approximation, reasonably good 
operation can always be obtained by locating the focus at the aperture 
of the horn. 
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CREEP ANALYSIS OF ROTATING SOLID DISKS WITH VARIABLE 
THICKNESS AND TEMPERATURE 


BY 


B. M. MA! 


SUMMARY 


A creep analysis for stress distributions in rotating solid disks of variable thick- 
ness and variable temperature is presented. This analysis is based on the theory of 
the Tresca criterion and the associated flow rule. The chief advantage of this 
analysis for stress distribution is that the complicated problem of disk design in gas 
turbines and jet engines can be solved with great simplicity by the closed form. The 
results obtained for the disk design will be on the conservative side. 

In this paper, the exponential function creep law at steady-state conditions will 
‘be employed. The profile of the meridian section of the disk is assumed to be of the 
hyperbolic form, and the temperature variation and the creep-rate distribution are as- 
sumed to be power functions of the radius of the disk. Some examples will be worked 
out for calculating the theoretical stress distributions in variable thickness and 
variable temperature disks. Stress distributions in uniform-thickness and variable- 
temperature disks will also be given for the purpose of comparison. 


NOMENCLATURE 


The following nomenclature is used in the paper: 


c, ¢, l, m, K = constant 
h = f(r) = cr™ = cro™x™ = hox™ = thickness of disk 
\ = a multiplier or a parameter which is a function of x 
r = radius to any point of disk 
r; = inner radius of disk 
= outer radius of disk 
= radius to the point x, of disk 
= time 
= absolute temperature 
= rate of radial displacement of disk at radius r 
= ri/To 
= re/To 
Xm = the point at which o, and o, are maximum 
=> 
= 
= 
= 
= 
= 


o = equivalent or effective stress 
oo. = pw ro = equivalent inertial stress in rotating thin ring of radius 7, 
oo = material creep modulus, having dimensions of stress 
1, Tr, 0, = tangential, radial and axial stresses of disk, respectively 
or, = radial stress at 7, of disk 


radial stress at x, 
tangential stress at x, 
maximum tangential and radial stresses, respectively 


Tiny 
creep rate or effective creep rate 


€ 


1 Department of Mechanical Engineering, South Dakota State College, Brookings, S. D.; 
currently on leave of absence and with the Department of Mechanical Engineering, University 
of Michigan, Ann Arbor, Mich. 
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€:, €r, €&¢ = tangential, radial and axial creep rates of disk, respectively 
€t, €r, € = tangential, radial and axial creep strains of disk, respectively 
p = y/g = density of disk material 
= specific weight of disk material 
g = gravitational acceleration 
w = angular speed of disk 


1, INTRODUCTION 


This paper may be considered an extension of the previous paper 
(1)? in which some formulas were derived for creep strains and stress 
analysis in rotating solid disks of variable thickness and uniform tem- 
perature. In the previous paper the temperature variation from the 
rim to the center of gas-turbine and jet-engine disks was not taken into 
consideration. Consequently, the mathematical derivations of the 
strain and stress formulas were less involved. 

Actually, the disk itself carries the turbine blades. Both the blade 
roots and the rim of the disk are highly heated by contact with hot 
gases and by conduction from the turbine blades. Consequently, 
higher temperatures and steeper thermal gradients are developed in 
these parts than those near the central portion of the disk. At the 
same time, the creep effects are cumulative near the periphery of the 
disk. Therefore, the combine factors of relatively high temperatures, 
steep thermal gradients and cumulative creep effects will usually cause 
the creep of the disk to be most severe where the blade roots are held. 

Since the rotating disks of gas turbines and jet engines are operated 
under high temperatures and high speeds, the centrifugal forces may 
produce stresses that exceed the yield strength of the disk material with 
the result that plastic flow is produced. 

The interaction of this plastic flow and creep, together with the 
variation of thermal gradients at high temperatures in the disk, may 
cause further serious stress distributions and even failures that might 
not be predicted on the basis of purely elastic- or plastic-stress analysis. 

In this paper, the stress analysis of rotating solid disks having 
variable thickness and variable temperature will be presented. This 
analysis will be based on the theory of the Tresca criterion and its 
associated flow rule, and will use the exponential function creep law 
at steady-state conditions. Examples which illustrate the application 
of the theoretical analysis for stress distributions in the disk will also 
be given. 

In general, the stress distributions in the disk may be divided into 
the following two regions :* 


(1) o, >o, >o, = 0,x. < x <1 for the outer region of the disk 


(2) o, =o, >o0, =0,x; < x < x, for the inner region of the disk 


2 The boldface numbers in parentheses refer to the references appended to this paper. 
3 The symbols used in this paper are defined in the Nomenclature Section on pp. 1-2. 
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It should be noted that in the inner region of the disk the pattern of 
stress distribution is identical to that in the disk of uniform strength, 
namely o, = o, = constant. 

A solid disk of uniform strength which may be called the exponential 
disk is suitable for relatively higher peripheral velocities, ¢,./o., (when 
the material and outside radius of the disk are given) of high-speed gas 
turbines and jet engines. It may be noted that the hyperbolic disk 
assumed in the present paper gives the same pattern of stress distribu- 
tion in the most-stressed region of the disk as the disk of uniform 
strength, but has the advantage of less mathematical involvement than 
the exponential disk. 


1.0 


(Cc, s 1.0) 


f,(x)=c,x! 


0.4 0.6 


X=r/to 


Fic. 1. Temperature distributions in rotating solid disk. 
2. THE BASIC EQUATIONS 


As discussed in a previous paper (2), the creep rate, é, is assumed to 
be equal to the product of a function of effective stress, F(a), a function 
of time, f(t), and a function of temperature, g(7), that is, 


é= f()g(T). (1) 


It is further assumed that the function of effective stress, F(¢), may 
take the form of the exponential funciion creep law, e*/*+, for relatively 
high stresses or high temperatures. Where g is the effective or equiva- 
lent stress in the maximum shear theory and ¢, the material creep 
modulus. As to the function of temperature, g(7), in taking the tem- 
perature variations into account, we further postulate that the tempera- 
ture varies with a power function of the radius 7 or x of the disk (3), 


3 
| 
0.6 HY 
0.2 Ly 
0.2 0.8 1.0 
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and the creep rate likewise increases with the power function of x at the 
steady-state condition, provided that F(c) remains unchanged. Hence 


= gff(x)] = filx) = (2) 
é= F(o) f(t) filx) (1a) 


where c, and / are constants. The function f,;(x) for the various values 
of / and x (c; = 1) is plotted in Fig. 1. 

Based on the foregoing discussion, the basic assumptions given in 
the previous paper (1), and the associated flow rule of Tresca (4), 
Eq. 1a for the problem of disks may be written as follows: 


and 


+ otlte 


é- = — =0 


< 1; and 


foro, > 0, >Oand x, < x 


ou = «= (4) 


— (& + é) = — (t) fi (x) 


for o, = 0, >0, x; < x < x. and $ <X < 1; where J is a multiplier 
or a parameter which is a function of x. Here the incompressibility 
equation, é, + é, + é€. = 0, has been used in Eqs. 3 and 4. 

Integrating Eqs. 3 and 4, the basic equations for creep strains are 
obtained. 


— = Kerf, (x) [soa (5) 
0 


€, 


forg, >¢, >0,x,. <x < 1l;and 


AKer f(x) f 


(1 — (x) [roa (6) 


2 
. 
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fora, = 0, >0,x; <x <x,and} <2 <1. The initial creep strains 
have been assumed to be zero. It should be also noted that at the 
steady-state condition the effective stress is independent of the time. 

The well-known equilibrium equation for the element at radius 
r or x of the rotating disk, Fig. 2, is given by 


d(hxo,)/dx — ho, + o,hx*? = 0 (7) 


Fir. 2. Rotating solid disk. 


in which / is the variable thickness of the disk and o, = pw*r,?, the 
equivalent inertial stress in rotating thin ring of radius r.. 

In order to determine the stresses o, and ¢,, it is necessary to intro- 
duce one more independent equation in addition to the equilibrium 
equation. This is the equation of compatibility since the disk is one 
of the redundant structures. 

Thus the known compatibility equation for creep rate of the disk 
may be simply given by 


ax 


3. DERIVATIONS FOR THE STRESS EQUATIONS OF CREEP 
1. Outer Region of the Disk 

Substituting Eqs. 3 into Eq. 8 and simplifying, we have 
, dfi 


1 


4 

6 6 
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dr 
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Upon integration between x, and x, this yields in the non-dimensional 
form 


_ Foy Xf i(%e) 


xfi(x) 


(9) 


in which o;,, is the tangential stress at x = x.. Now introducing Eq. 9 
into Eq. 7 and integrating between x, and x, the expression for the 
radial stress is obtained. 


where ¢o,, = at the point x = x,. Nowif we let x = l anda, = a,, 
at the outside radius of the disk r,, and then substitute these values into 
Eq. 10 and rearrange, the following equation for stress at x, is found. 


(ia. +f hdc). (11) 


Upon substitution of the values h = cr,"x" = h,x™ for the hyperbolic 
profile of the disk (c and m are constant but m is a negative value), and 
fi(x) = cix', fi(x.) = cx,' for the temperature or creep-rate distribu- 
tions in the disk. Equations 9, 10 and 11 for x, < x < 1 may become, 
respectively, 


Ca Ta x 


1 
1+m Ca 


1+m 1+m Ca 3+m (10a) 


—y ltm — 3+m 


It can be readily verified that at x = x., Eqs. 9a and 10a reduce identi- 
cally to o:, = ¢,,, which can be determined from Eq. 11a when the 
values of 1, m, x., ¢./¢. and a@,,/o, are prescribed. It is also noted that 
for the particular case when x, = 1, Eq. 11a reduces to o,, = o,, = Or» 
at 
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Now if we assume values of / = 1.5 and 2.5, m = — 3, r, = 10 in. 
and ¢,/o, = 0.0275 (this corresponds to the case where ¢, = 7,000 psi., 

= 254,500 psi. at a disk speed of 17,800 rpm., and y = 0.283 lb. per 
cu. in.) and substitute these values in Eq. 11a, the resulting values of 
¢,,/0, plotted against x. for various values of ¢,,/o_ will be as shown 
in Figs. 3 and 4. 

Either the inertial radial stress of the disk, ¢,,, or the total centrifugal 
force of the turbine buckets or blades, that is transmitted through the 
rim of the disk to the disk proper, can be approximately calculated for 
any given dimensions of the disk and blades. To calculate this, we 
may integrate the elemental centrifugal forces from the tip to the root 
of the blades. In addition, the centrifugal force due to the rim of the 
disk may be also included. 

In order to determine the critical value of x (=x.) at which the 
condition ¢, > > 0 will not be satisfied, we divide Eq. 10a by 
Eq. 9a; differentiate the ratio o,/o, with respect to x and put it equal 
to zero for a maximum at x = x,. This results in 


/ 
£5 = | (taleg the positive eign): (12) 


With the given values of 1, m, ¢./a. and o;,/¢, as calculated from 
Eq. 11a or selected from Fig. 3 or 4, the critical value of x, can be readily 
found. Take for example, / = 1.5,m = — 3}, ¢,,/o. = 0.4 (the largest 
value of ¢,,/¢, to be assumed in the later examples) and o;,/¢. = 0.568 
from Fig. 3. The resulting critical value of x. as calculated from Eq. 12 
(using the above values of /, m, o:,/¢, and @,./¢, = 0.0275) is 0.595. In 
this particular case, any assumed value of x. less than the calculated 
value of 0.595 will not satisfy the required condition o, > o, > 0, over 
a portion of the outer region of the disk. This is the reason why the 
value of x, will be assumed to be 0.60 or larger in the curves shown in 


Figs. 5 and 6. 


2. Inner Region of the Disk 


In this region, since o, = o,, Eq. 7 can be simply integrated between 
x and x,., we have 


Ga Ga 


Introducing the thickness function 4 = h.x" and h, = h.x.", now we 
have 


Ca Oa ( x ) Oa (2 + m)x™ (13e) 
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To find the maximum or peak stresses in the disk, which is of most 
interest in practical design, we differentiate Eq. 13a with respect to x, 
set it equal to zero and solve for x = x, where the maximum stresses 
occur. This gives 


Xn = ig, "+5 (14) 
t (2 + m)xn™ (15) 


In order to determine the value of X if it satisfies the condition 
4} <d <1 forx,; < x < x., we now substitute Eqs. 4 into Eq. 8, differ- 
entiate and rearrange. Thus 


Using f: = fi(x) = cx', integrating this equation between x and x., 
and taking \ = 1 at x = x,, we obtain 


x, \ 2+! 1 ty 
Az={— — er tidy (16) 


x 


in which o;, is given by Eq. lla and « is given by Eq. 13a. By a 
numerical evaluation of the integral in Eq. 16, it may be found whether 
or not the condition } <A <1 holds for x; <x <-x,.. Experience 
shows that this is usually true. 

The procedure for stress analysis in rotating solid disks of variable 
thickness and variable temperature may be summed up as follows: 


1. For the outer region of the disk, ¢, > ¢, > 0 forx, <x < 1. 


(a) Equation 11a is used to determine o,,/o, for given values of /, 
Mm, X-, Go/¢_ and o,,/a_ (as shown in Figs. 3 and 4). The value of o, 
is obtained from tension creep test data of the disk material at elevated 
temperatures, and ¢,, is obtained from the total calculated centrifugal 
force of the turbine blades transmitted through the rim of the disk 
to the disk proper. 

(b) The value of o;,/0. as determined in step (a) is next substi- 
tuted into Eq. 12 and the critical value of x, is found and checked 
against the assumed value of x, in Eq. 1la (as shown in Figs. 3 and 4 
or similar charts). 

(z) If the value of x, found by Eq. 12 is smaller than or equal 
to that assumed in Eq. 11a, then the condition o, > ¢, > Oforx, <x <1 
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is satisfied and this value of ¢,,/¢, and the assumed value of x, in Eq. 1la 
can be used for further calculations or follows. 

(it) If the value of x. found by Eq. 12 is greater than that as- 
sumed in Eq. 11a, then the above condition will be not satisfied and 
steps (a) and (b) must be repeated, assuming a larger value of x, in 
Eq. 11a, until the value of x. found by Eq. 12 is smaller than or equal 
to that assumed in Eq. 11a. In this particular case, it would be con- 
venient to have Figs. 3 and 4 or similar charts for guidance, to determine 
o,,/¢_ and the corresponding assumed value of x.. 


0.7 
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Fic. 7. Creep stresses in rotating solid disk, / = 2.5, m = — 0.5, x. = 0.7. 


(¢) Introducing these values of x, and o:,/o. from (6) to Eqs. 9a 
and 10a, the stress distribution of ¢,/¢, and ¢,/o, are then determined 
and the condition o, > ¢, > 0 for x. < x < 1 is satisfied. 
<x < x, and 


2. For the inner region of the disk, o, = o, > O for x; 
$<A <1. 


(a) Equation 13a is used to compute ¢,/o. = ¢,/o, after the values 
of o,,/¢, and x, are substituted from step 1b above. 
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(6) The value of \ is now checked by Eq. 16 to determine if the 
condition } <\ < 1 for x; < x < x, holds. 

(c} Finally, the value of x, is calculated from Eq. 14 and the 
maximum stress, ¢:,,/¢e = Orp,/Ge at Xm are determined from Eq. 15. 


From the procedure described above, the stress distributions in the 
variable-thickness and variable-temperature disks having / = 1.5, 2.5, 
m = — 3 and o./o, = 0.0275 (at temperatures near 1200° F.) are cal- 
culated and plotted against x for the various values of x. and ¢@,,/aa. 
The results are shown in Figs. 5, 6, and 7, respectively. 

For purposes of comparison, the stress distributions of the uniform- 
thickness (m = 0) and variable-temperature disk (other data remain 
the same) are plotted against x in Figs. 8, 9, and 10, respectively. 
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Fic. 10. Creep stresses in rotating disk, 1] = 2.5, m = 0, x. = 0.7. 
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Furthermore, the stress distribution of the variable-thickness and uni- 
form temperature disk at x. = 0.6 is given in Fig. 11, and the elastic 
stress distribution of the variable-thickness disk without thermal effect 
is given in Fig. 12. In calculating the elastic stresses of the disk the 


following equations are used : 


__ 8+(3+)m -91—(1+m) _1+3n | 


in which m= —3; g,=0.851, one of the roots for the hyperbolic disk ; 
nu = 0.3, Poisson’s ratio and ¢,,/¢, = 0, 0.2 and 0.4, respectively. It 
is interesting to note that for this special case, the peak radial stresses are 
invariably greater than the peak tangential stresses in the elastic disk. 

The -values of x,, and (=@,,,/@a) Calculated from Eqs. 14 and 
15 for Figs. 5, 6 and 7 are given in Table I. 


TABLE I.—Maximum Stresses o1,,/0a = in the Rotating Solid 
Disk of Variable Thickness and Variable Temperature. 


m = — 0.5 
Xe l Ses = 
o:,/¢. 0 0.2 0.4 0 0.2 0.4 
0.6 1.5 0.3791 0.4595 0.5351 0.2873 0.4238 0.5726 
0.3851 0.4640 0.5367 0.2963 0.4305 0.5761 
2.5 0.4645 0.5370 0.2976 0.4316 0.5769 


By comparing Figs. 5, 6 and 7 with Figs. 8, 9 and 10, we see that the 
maximum stresses produced in the variable-thickness and variable- 
temperature disk are much less than those produced in the uniform- 
thickness and variable-temperature disk as we would expect. 


4. DISCUSSION AND CONCLUSIONS 


From the foregoing creep analysis, the following discussion and 
conclusions may be made. 

1. It is seen from Figs. 5, 6, 7 and Table I that the maximum stresses 
in the variable-thickness and variable-temperature disk increase with 
the centrifugal force or inertial stress, ¢,,/¢2, of turbine blades, and the 
peaks of the maximum stresses have a tendency to shift outward from 
the center of the disk as the values of o,,/o4 increase. 

2. Temperature variations or thermal gradients tend to increase the 
maximum stresses moderately (comparing Figs. 6 and 11). The maxi- 
mum stresses, however, appear to be insensitive to the changes of the 
temperature variation (from / = 1.5 to / = 2.5 in this particular case) 
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as the turbine speeds, or the values of ¢,,/o. increase (Figs. 6 and 7). 
The thermal effects on the maximum stresses gradually diminish and 
probably vanish as o,,/o4 increases to the burst speed of the disks. 

3. The maximum stresses increase very slightly with x,., Figs. 5 and 6 
or Table I. Therefore, the value of x. found in Eq. 12 either smaller 
than or equal to that assumed in Eq. 11a only has slight effect on the 
maximum stresses. 

4. There are two possible extreme cases of stress distributions for 
variable-thickness and variable-temperature disks. This also applies 
to uniform-thickness and uniform-temperature disks and to variable- 
thickness and uniform-temperature disks (1). If 

(a) o, >o, > 0 at x; < x <1, then only Eqs. 9a and 10a are 
required to determine the stresses except at the point of the disk center. 

(6) o, =o, >0 at x; < x <1, then only Eq. 13 is needed to 
calculate the stresses, and the stress pattern is exactly like that in the 
disk of uniform strength (this may be called the ‘‘pseudo-exponential 
disk’). 

5. It is seen from this creep analysis which is based on the Tresca 
criterian and its associated flow rule, that by using the exponential 
function creep law, the incompressibility equation, the equilibrium 
equation and the compatibility equation; the stress equations were 
derived with simplicity and ease. In addition, this stress analysis 


appears to be on the conservative side for the design of gas-turbine and 
jet-entinge disks, since the effective stress in the Tresca criterion is 
higher than that in the Mises criterion. 
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A TECHNIQUE FOR STUDYING SUPERSONIC COMBUSTION 


A technique for securing information on the combustion process at 
supersonic speeds has been developed by the National Bureau of Stand- 
ards with the support of the Navy Bureau of Weapons and the Air 
Force Office of Scientific Research. The technique involves firing a 
missile at hypersonic velocities through a combustible gas. High-speed 
schlieren photography is used to obtain data concerning shape and 
detachment distance of the shock wave, gas density, and other effects. 
These data are evaluated to derive information on the combustion 
process in the vicinity of the missile. 

This experimental method, devised by F. W. Ruegg of the combus- 
tion controls group, permits precisely controlled variation of such param- 
eters as pressure, temperature, fuel, and speed over a wide range under 
realistic conditions and at very little cost. It can be employed to ex- 
amine in advance some of the problems of stabilizing detonation waves 
on a stationary object in a moving combustible mixture. If the wave 
can be established at the front of the missile, the drag will be substan- 
tially reduced. Results of studies using this method might be applied 
to the generation of thrust by burning at supersonic velocities either 
on the exterior of a missile or within it. Also, the method should con- 
tribute to the design of flame holders for supersonic combustion. 

During the investigation, laminar combustion in a supersonic stream 
was observed for what is believed to be the first time. An unusual 
occurrence was that of large-amplitude oscillations in the megacycle 
frequency range produced by the combustion process at the front of 
the missile. 

The study of the combustion process in standing detonation waves 
and in supersonic streams is of great practical importance in the field 
propulsion of hypersonic missiles and planes. Extensive information 
obtained in tube experiments is available concerning moving, plane 
detonation waves, but this body of data cannot with any assurance be 
simply transposed to the converse situation of a moving gas and a 
stationary oblique wave on a solid object. Furthermore, the object 
itself may exert unforeseen influence. 

Some information has recently been obtained on the action of both 
normal and oblique shocks in inducing combustion in closed ducts and 
in free jets. The results are somewhat ambiguous, however, because 
of fuel and air mixing problems that exist in streams at high speed and 
temperature. Also, the range of experimental conditions was limited 
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by difficulties encountered in attempting to produce high-speed streams 
at high stagnation temperature and pressure. 

The technique developed at the Bureau overcomes these obstacles by 
creating conditions such that a shock wave is established ahead of a 
high-speed missile, and then determining the conditions required to in- 
duce combustion and observing the effects of the missile on the com- 
bustion process. The stagnation temperatures attained are realistic, 
and the results are not affected by the method of mixing air and fuel 
or by turbulence. The experiment can easily be modified to determine 
the effects of the temperature and pressure history of the combustible 
gas on the combustion in the shock wave. 

The Bureau’s experiments were performed at Mach numbers from 
4.2 to 6.5 in a mixture of hydrogen (30 per cent by volume) in air. 
Combustion was initiated by means of compression and heating within 
the shock wave standing on the missile. Characteristics of the combus- 
tion process depended strongly on the Mach number and on the state 
of the gas ahead of the wave. 

The experimental apparatus comprises three essential parts: a gun 
and a range that contains the missile and the test gas, an electronic 
time-delay-ratio generator that computes the time to fire a high-voltage 
flash lamp, and an optical system for schlieren photography. 

In the experiments performed up to the present, 0.785-inch-diameter 
nylon spheres were used as missiles. The sphere is mounted in a 20- 
millimeter case containing a fast-burning gun powder and fired from a 
.60-caliber, 6-ft.-long test barrel through a blast chamber and the test 
section to a heavy metal target and trap. Two phototubes are located 
upstream of the test gas, which is isolated from the remainder of the 
system by two plastic diaphragms. 

Elapsed time of flight of the missile between the two phototubes is 
measured by an electronic timer that counts the number of pulses from 
an oscillator operating at a frequency of one megacycle. Signals from 
the phototubes enter the time-delay-ratio generator, which computes 
the time of arrival of the missile at the photographic station and triggers 
the high-voltage flash lamp at that time. A second electronic timer 
measures the time between the missile’s arrival at the second phototube 
and the time at which the lamp is flashed. Two speed measurements 
are thus available for calculating the deceleration of the missile in the 
test section. 

HIGH-RESOLUTION CAMERA 
A research camera capable of projecting a parallel line pattern of 
50,000 lines per inch has been developed at the National Bureau of 
Standards. The camera, designed by C. S. McCamy of the Bureau’s 
photographic research laboratory, is a step toward the development of 
a standard method for determining the resolving power of photographic 
materials. 
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At the present time, there is no standard method for measuring 
resolution. The Bureau's work is part of efforts by a committee of the 
American Standards Association to establish such a standard. Five 
laboratories in the United States and Canada are working on methods 
which can be employed. By comparison of results, they will be able 
to determine the factors which must be specified in order that all reso- 
lution tests on the same film or plate may yield the same result. 

To insure that any limitations revealed by the test are those of the 
photographic materials and not of the test equipment, the camera must 
be able to produce a line pattern much finer than can be resolved by the 
emulsion being tested. To test the new camera’s full capabilities 
(50,000 lines per inch) it was necessary to use extremely high-resolution 
spectrographic plates having a film speed of 0.025 when rated on the 
ASA scale for pictorial negative materials. 

The camera’s optical system is essentially a fine apochromatic micro- 
scope operated in reverse to make an extremely small image of a large 
chart. The camera has interchangeable lenses. Objectives in use are 
an 8-mm. lens having a numerical aperture of 0.65, which is equivalent 
to an f-number of 0.77, and a 16-mm. lens having a numerical aperture 
of 0.30 or an f-number of 1.66. For this application the elements of 
the objectives were respaced and the interiors of the lens mounts were 
specially blackened to minimize reflections. 

Film or plates to be tested are placed against a surface attached to 
the objective lens mount by a fine screw thread, with a gear mechanism 
permitting minute adjustment of focus. Films are held in place by a 
miniature vacuum back. The focus is first set by direct visual observa- 
tion with a microscope. Several exposures are then made, setting a 
slightly different focus each time, until the position of best focus is 
found. The camera is operated at a constant temperature to avoid 
changes in focus resulting from thermal expansion or contraction. 

To secure even illumination of the targets to be photographed, a 
compensator plate is placed between the lamp and condensers of the 
camera’s illuminating system. The compensator is a glass sheet with a 
circular pattern printed on it—heavier at the center than at the edges 
—which counteracts the tendency of the condensers to give higher 
illumination at the center of the beam. Provision has been made to 
place filters, either color temperature-correcting filters or neutral density 
filters, between the illumination and the test chart. Color correcting 
filters are used to adjust the illumination for the particular emulsion 
being tested. Neutral density filters are used to obtain a series of 
exposures of varying illuminance. The maximum resolving power ob- 
tained in the series is reported. 

One important consideration in photographing the fine line patterns 
is to eliminate mechanical vibrations which could blur the photographs. 
To this end, the camera, including test chart, illuminator, shutter, and 
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optics is encased in a massive steel cylinder. The cylinder, in turn, 
is suspended from the ceiling by springs. The vibrations naturally 
present in the building were determined to be mainly 30 cps, 60 cps, 
and 120 cps; since the springs used for support have a natural frequency 
of approximately 2 cps, an effective mechanical insulation is achieved. 

The new camera was used in an experiment to test the generally 
accepted correlation between resolving power and legibility at an ex- 
treme reduction ratio. In addition, the photograph made shows the 
extremely fine detail which may be recorded by this system. 

In this experiment, Bernard Fouquet of the Bureau’s photographic 
research laboratory made a photocopy of the first page of the Bible at 
a linear reduction ratio of 1000, thus reducing the area by a total factor 
of one million. 

The page was first reduced 2.5 times to fit into the camera’s chart 

holder; then it was reduced by a factor of 400 with the new camera. 
If the entire Bible could be reproduced at this reduction ratio, it would 
cover no more than } cm?, less than the area of Lincoln’s head on a 
penny. 
The microcopy was made on a high-resolution spectrographic plate 
using an exposure of 10 sec at f/0.77. It is clearly legible when viewed 
directly at a magnification of 1250 through a fine optical microscope. 
Blue light is used to view the microcopy since a short wavelength en- 
hances resolution. 

The Bureau does not recommend this method for routine micro- 
copying. Because of the difficulty in focusing, the problem of keeping 
dust from the field, the long exposures involved, and the difficulty of 
viewing the copy, the method is too tedious for routine use. 

Since one of the main purposes of the undertaking was to develop a 
standard method of measuring resolving power of photographic ma- 
terials, the materials and components used were deliberately limited 
to those generally available. Somewhat better performance might 
be realized if this limitation were relaxed. Much more work remains 
to be done before a method can be standardized. 
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STATISTICAL THEORY OF COMMUNICATION, by 
Y. W. Lee. 509 pages, diagrams, 6 X 9 in. 
New York, John Wiley & Sons, Inc., 1960. 
Price, $16.75. 


This book, written primarily for the first- 
year graduate student in electrical engineering, 
presents clearly and rigorously a physically 
motivated and systematic account of the 
statistical theory of communication. The 
author has pioneered in the synthesis of 
electric networks by means of the Fourier 
transforms of Laguerre’s functions and is 
co-inventor (with Dr. Norbert Wiener) of 
such networks in the time domain and the 
frequency domain. Dr. Lee and Dr. Wiener 
set forth for the first time in the history of 
circuit theory the amplitude-phase relations 
in system functions, and the real part and 
imaginary part relations in these functions, 
in the form of Hilbert transforms. So, when 
the author says in his Preface that he does 
not attempt to cover the subject of study but 
attempts to uncover it for the student, there 
comes a stimulating book on the fundamental 
theory of communication. 

The reviewer was aware of the intimate 
cooperation between Wiener and Lee when 
he was able to invite Dr. Wiener to visit 
Tsing Hua University for the academic year 
1935-36 so that they could work together. 
In his book “J Ama Mathematician” (Double- 
day and Company, Inc., Garden City, N. Y., 
1956), Wiener said: “If I were to take my 
specific boundary point in my career as a 
journeyman in science and as in some degree 
an independent master of the craft, I should 
pick out 1935, the year of my China trip, as 
that point.”” After World War II, Dr. Lee 
was invited to teach at M.I.T., and he has 
stayed on teaching and doing research since 
1946. Quoting Wiener again: ‘“‘What made 
Lee's situation easier was that I had recently 
developed a considerable part of the statistical 
theory of communication engineering .. . . 
Thus, he took over the problem of working out 
in detail the communication and engineering 
consequences of ideas which I had only 
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sketched in general terms and of making him- 
self interpreter to the engineering public of 
the field which I was later to call cybernetics.”’ 

The book starts with a short introduction 
which refers to the classic work of Wiener 
entitled ‘‘The Extrapolation, Interpolation, 
and Smoothing of Stationary Time Series with 
Engineering Applications’ (John Wiley & 
Sons, New York, 1949). Chapters 2 to 8 
give the basic concepts and tools which are 
summarized in Chapter 9. These chapters 
cover generalized harmonic analysis, proba- 
bility, random variables, ensembles, and dis- 
tributions, time averages and ensemble aver- 
ages, correlation functions and power density 
spectrums of random processes. In Chapter 
9, the author said: 


Messages and noise are not completely 
predictable so that unique specification 


of these functions for all values of time 
is contrary to their nature.... 
meaningful and useful method of descrip- 
tion has been found to be the description 
by a set of probability densities (proba- 
bility distributions if the ensemble is 
discrete). This is a fundamentally dif- 
ferent concept in describing functions as 
compared with the concept of complete 
determinism and the concept of ampli- 
tude and phase spectrums for periodic 
and aperiodic functions. 


The next eight chapters cover such impor- 
tant topics as sampling theory, detection of 
a periodic signal in noise, optimum linear 
systems, optimum filtering and prediction, 
errors in optimum systems, and theory of 
optimum systems. The fundamental rela- 
tions for a linear system discussed in Chapter 
13 include the input-output relation in auto- 
correlation, the input-output crosscorrelation 
theorem, and the crosscorrelation of linearly 
transformed random functions. In Chapter 
17 optimum systems are studied both on 
the minimum-integral-square-crosscorrelation- 
error criterion and with transient inputs on 
the minimum-integral-square-error criterion. 

The last two chapters—Chapters 18 and 
19—deal with the application of orthonormal 
functions to the representation of correlation 
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functions and power density spectrums and 

to the synthesis of optimum linear systems. 

As the author said in the Preface: 
. . . However, because of the great im- 
portance of the orthonormal representa- 
tion of a linear system in the statistical 
theory of optimum nonlinear systems, as 
discussed by Professor Wiener in his 
Nonlinear Problems in Random Theory, 
the last chapter of the book is devoted to 
the synthesis of linear systems by ortho- 
normal functions. 


When the reviewer was Visiting Professor 
at M.I.T. in 1950-52, Professor Lee was 
already well established in this new field and 
had started to write this book. It has taken 
him a rather long time to finish the manu- 
script as he wanted to present the material 
as clearly and rigorously as possible. It is 
the reviewer's great privilege to recommend 
this carefully written treatise on the modern 
statistical theory of communication to the 
engineering public. The author has ably 
demonstrated the importance of this new 
field which is ‘a landmark in the history of 
communication theory.” 

Y. H. Ku 
University of Pennsylvania 


Jet PropuLtsion ENGINEs, edited by O. E. 
Lancaster. 799 pages, plates 6 X 9} in. 
Princeton, Princeton University Press, 
1959. Price, $20.00. 


This book consists of a group of well inte- 
grated articles on jet propulsion. Beginning 
with its historical development, chapters on 
general formulas, cycle analysis, matching 
components, engine performance, with inter- 
esting chapters on fuels, ram rockets, heli- 
copter jet rotors and atomic energy applica- 
tions cover a wide range of technical informa- 
tion on this subject. 

The entrance chapter by F. J. Malina gives 
an interesting discussion of the history of 
developments, including R. H. Goddard's 
pioneer experiments. It is impressive to note 
his introduction of the converging-diverging 
nozzle for increasing jet speeds above the 
acoustic velocity, enormously upping the ef- 
ficiency, with ejection velocities increased 
from 1000 to over 8000 ft/sec. 

Chapter I by Maurice Roy covers a detail 
analysis on the interaction of flows with the 
engine and airplane exterior. The effective 
thrust is divided correctly into interior and 
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external components. Next the momentum 
change of the flow through engine is de- 
termined, where the housing reaction and end 
pressures at exit and entrance are considered. 
The latter are determined by an analysis of 
entrance and exit flow conditions. These are 
complicated by the interaction of internal and 
external flow which takes place at the ends, 
as well as the appraising of the latter, and its 
drag reaction on the exterior of the engine 
housing. 

In a succeeding Chapter C, C. A. Meyer 
follows a similar approach, along with discus- 
sions of experimental data on axial and cen- 
trifugal compressors, burners, and perform- 
ance matching. The experimental data and 
discussions are worthy of study. Section D 
by Driggs and Lancaster covers very useful 
studies relating to performance analysis, pro- 
peller performance, turbine design, controls, 
gear box and accessories, as well as engine 
testing methods. 

The general discussion on ram jets by 
DeMarquis, D. Wyatt and Bruce T. Lundin 
is extensive and most interesting. Similar 
studies are not likely to be found elsewhere. 
Special chapters on intermittent jets by 
Joseph V. Foa include propulsive cycles, the 
pulse jet and wave engines. 

Section G by Martin Summerfield on the 
liquid propellant rocket engine is a valuable 
discussion of basic thermodynamic relations 
as to the nature of combustion and nozzle 
performance with liquid fuels. It includes 
cooling problems of rocket motors. 

Section H, by C. E. Bartley and Mark M. 
Mills covers the interior ballistics of solid 
propellants. This important subject is well 
treated and gives features on mechanical 
design as well. 

The Ram Rocket, by Irvin Glassman and 
Joseph V. Charyk, is of special interest in its 
relation with the ram jet, as well as the general 
problem of mixing and combustion of two or 
more dissimilar streams. The work covers 
considerable analysis as well as experimental 
performance evaluation. 

The chapters by Albert Gail on Jet Rotors 
cover the kinematics and blade propulsion 
reactions with jet motors. This extends to 
helicopter jet rotor performance, including 
mechanical design. 

A final chapter on the application of atomic 
energy in jet propulsion is covered in a com- 


2 
4 
gg 
: 


62 Book REVIEWS 


prehensive way by Ralph Zirkind. Reactor 
design, shielding and general systems of heat 
transfer through liquid metals to turbo jet 
radiators are accompanied by preliminary 
design data. 

A closing chapter by F. Zwicky discusses 
the future prospects of jet propuslion. 

In spite of the outstanding array of very 
worthwhile technical discussions, covering a 
fair range of the subject, it would seem an 
over-all basic analysis on the reaction of 
fluids, particularly as to the nature of energy 
transfer and thermodynamic relations between 
the flow and engine housing and jet, and the 
significance of propulsion efficiency, should be 
included in a condensed entrance chapter. 

On the whole, this book should prove of 
great value to all technical engineers con- 
cerned with jet propulsion. 

RUPEN EKSERGIAN 
The Franklin Institute Laboratories 


REACTOR ANALYsISs, by Robert V. Meghreb- 
lian and David K. Holmes. 808 pages, 
illustrations, 6 X 9in. New York, McGraw- 
Hill Book Co., Inc., 1960. Price, $19.50. 


This volume is by far the most complete and 
extensive treatise on the analytical methods 
used in the study of the properties and be- 
havior of nuclear fission reactors. The con- 
tent of the book is an extension of the notes 
developed by the authors for a course given 
by them for several years at the Oak Ridge 
School of Reactor Technology (ORSORT) at 
the Atomic Energy Commission's Oak Ridge 
National Laboratory. 

The subject of reactor analysis deals with 
the methods and mathematical models used 
in predicting the spatial, temporal and energy 
distributions of the neutron population in 
assemblies containing fissionable substances. 
Since the authors emphasize the methods of 
analysis rather than the physical processes 
themselves, this test serves to complement 
another recent book which places its major 
emphasis on the physical aspects of fission 
reactors (‘‘The Physical Theory of Neutron 
Chain Reactors,” by A. M. Weinberg and 
E. P. Wigner). It should be noted that both 
these books cover much the same subjects, 
and both present the descriptions of the phys- 
ical phenomena, the theories, and the formula- 
tions of methods of analysis; it is in their 
placing of emphasis in which they differ. 


Bik 


The authors have proceeded logically 
through the principal topics of reactor analy- 
sis. Beginning with an introductory dis- 
cussion which describes the various types of 
nuclear reactors and the various problems 
confronting the reactor physicist (Chapter 1), 
they proceed through extensive presentations 
of probability concepts (including nuclear 
cross sections) and neutron multiplication 
and flux (Chapters 2 and 3, respectively). 
There then follows a group of chapters with 
each devoted to a separate major topic: 
slowing-down in an infinite medium, diffusion 
theory, Fermi age theory, transport theory, 
reflected reactors, reactor kinetics, and hetero- 
geneous reactors. The book concludes with 
three chapters covering relatively advanced 
subjects: control-rods, hydrogenous systems, 
and perturbation theory. 

The discussion of each method or model is 
presented in considerable detail in order to 
emphasize the approximation and limitations 
attendant to each. The authors have at- 
tempted to present the most suitable applica- 
tions for each of the theoretical models dis- 
cussed, and have in a majority of instances 
demonstrated the general procedure by means 
of a numerical example. The subjects dis- 
cussed are well documented by means of 
footnotes. A group of problems associated 
with the various topics is also included fol- 
lowing many of the chapters. 

The authors have intended this text for the 
graduate engineering student. Although it is 
well suited as both an introductory and an 
advanced text on reactor analysis, a basic 
nuclear physics background is presupposed, 
and the student or reader acquainted with 
the various types of reactor physics preblems 
will better appreciate the utility of the dif- 
ferent analytical models and the circumstances 
under which each is most useful. The authors 
have succeeded in writing an exceptionally 
clearly written text which presents the theories 
and the basic methods available for deter- 
mining neutron-density distributions in var- 
ious reactor configurations. The book is 
best suited for those intending to specialize 
in reactor physics or reactor analysis, although 
anyone who has the proper background in 
engineering mathematics could profit from 
its study. 

W. H. STEIGELMAN 
The Franklin Institute Laboratories 
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THE Major ACHIEVEMENTS OF SCIENCE, by 
A. E. E. McKenzie. Vol. I, 368 pages, 
plates, diagrams, 54 X 84in. Price, $5.50. 
Vol. I], 195 pages, illustrations, 5} X 84 
in. Price, $3.50. New York, Cambridge 
University Press, 1960. 


This two-volume work covers the history 
of the main developments in science from 
ancient times to the present. In the 24 
chapters comprising Vol. I, the author de- 
scribes such general concepts as the Copernican 
theory, the circulation of the blood, the crea- 
tion of modern chemistry, the atomic theory, 
the wave theory of light, evolution, the clash 
between science and religion, the theory of 
relativity, the quantum theory and the latest 
theories about the creation of the universe. 
In addition to relating events leading up to 
the specific development, Mr. McKenzie gives 
background information on the discoverer of 
the development and describes the effect of 
the discovery on contemporary society. The 
author has skillfully blended his capsule views 
of the important discoveries and develop- 
ments in science into an over-all picture of 
their influence on Western thought and civili- 
zation. In Vol. II, Mr. McKenzie has col- 
lected 91 extracts from original writings of the 
scientists whose theories he expounded in Vol. 
I. Arranged under the same chapter head- 
ings as Vol. I, these original writings have 
been carefully chosen from books, lectures 
and letters rather than the highly scientific 
papers of the men, so that the general reader 
may obtain an insight into the personalities 
of the men responsible for our scientific 
advances. 


INFRARED METHODS: PRINCIPLES AND APPLI- 
cations, by G. K. T. Conn and D. G. 
Avery. 203 pages, diagrams, 6 X 9 in. 
New York, Academic Press, 1960. Price, 
$6.80. 


The authors have prepared this introduction 


to infrared studies primarily for the research 
student in physics, chemistry, boichemistry, 
biology and medicine, because 


- we 
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believe that there is no other working man- 
ual."" Dividing the book into two parts, the 
authors cover principles in Section A and 
practical applications in Section B. Part A 
includes sources of radiation, optical ma- 
terials, detectors, amplifiers and dispersive 
systems. Part B is devoted to calibration of 
detectors, simple monochromators, instru- 
ments for gas analysis, and radiation pyrom- 
etry. The book is restricted to the ‘near 
infrared’’—the range from 1 to 254. One 
of its best features is the inclusion of over 
500 references, divided among the chapters 
and covering the period from 1863 through 
1959. 


QUALITATIVE THEORY OF DIFFERENTIAL 
Equations, by V. V. Nemytskii and V. V. 
Stepanov. 523 pages, 6 X 9 in. Prince- 
ton, Princeton University Press, 1960. 
Price, $12.50. 


This translation, based on the enlarged and 
revised second Russian edition, is No. 22 of 
the Princeton Mathematical Series. Part 
One (306 pages), dealing with classical dif- 
ferential equations, has four chapters: Exis- 
tence and Continuity Theorems; Integral 
Curves of a System of Two Differential Equa- 
tions; Systems of » Differential Equations; 
and A Study of Neighborhoods of Singular 
Points and of Periodic Solutions of Systems 
of m Differential Equations. In addition, 
there is a translation of a resume written by 
Nemickii of work done on qualitative theory 
(not included in the Russian work). Part 
Two (217 pages) covers topological dynamics 
and ergodic theory. 


MopDERN NUCLEAR TECHNOLOGY, edited by 
Mark M. Mills, Arthur T. Biehl and Robert 
Mainhardt. 336 pages, illustrations, 6 X 9 
in. New York, McGraw-Hill Book Co., 
Inc., 1960. Price, $9.50. 


This book is comprised of a series of lectures 
given under the Engineering Extension Pro- 
gram of the University of California during 
the summers of 1956 and 1957. Approxi- 
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mately twenty-two guest speakers, all special- 
ists in the field, participated in the course 
entitled, “A Survey of Nucelar Engineering 
for Management.” The purpose of this sur- 
vey course was to present on a management 
level the fundamentals of nuclear engineering, 
in order to help management personnel under- 
stand some of its problems. Basic engineer- 
ing principles or method of approach for the 
solution of engineering problems are therefore 
not discussed. The presentation of material 
is then directed toward the mature industrial- 
ist with an engineering degree, but inactively 
engaged in engineering for five to ten years. 


ELEctTrRIcAL NoIsE, by W. R. Bennett. 280 
pages, diagrams, 6 X9 in. New York, 
McGraw-Hill Book Co., Inc., 1960. Price, 
$10.00. 


This book presents a comprehensive ac- 
count of how noise originates in electrical 
circuits, how it is measured, and how circuits 
may be designed to minimize undesirable 
effects from noise. It describes in qualitative 
terms the physical nature of various important 
noise sources, including thermal agitation or 
resistance noise, shot noise in vacuum tubes 
and semi-conductor junctions, noise from 
spontaneous emission of electromagnetic radi- 
ation, and noise in gas discharges. The funda- 
mentals needed for analyzing basic sources 
of noise are stressed as well as the methods 
of measurement and design. In addition to 
the standard theory of noise figure and its 
significance, a treatment of the more compre- 
hensive Haus-Adler theory of noise measure 
is given. The author also reviews noise in 
the various methods of signal transmission 
such as amplitude modulation, frequency 
modulation, and the different kinds of pulse 
modulation. 


MANAGEMENT CONTROL SysTEMs, edited by 
Donald G. Malcolm and Alan J. Rowe. 
375 pages, diagrams, 6 X 9in. New York, 
John Wiley & Sons, Inc., 1960. Price, 
$7.25. 


These are the proceedings of a symposium 
held at the System Development Corporation 
in Santa Monica, California, July 29-31, 1959. 
A comprehensive treatment on management 
controls by some thirty experienced individ- 
uals in all areas of management is presented 
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in this volume. The editors and contributors 
have constructed a case for “designed” con- 
trols to produce better results in management 
patterns now and for the long-range future. 
The problem of designing adequate controls 
is herein given careful consideration, especially 
in relation to the fact that the design of 
systems is currently at best a systematic 
process rather than one where analytical 
evaluation is possible. Special attention is 
given to computers and their role in this 
problem. Contributors to this symposium 
have described and weighed the advantages 
of using computer applications to simulate 
company activities so as to create a manage- 
ment-laboratory where the effects of policy 
and procedures can be tested and evaluated 
prior to their adoption. 


ComPuTER Locic: THE FUNCTIONAL DESIGN 
or DicitaL Computers, by Ivan Flores. 
458 pages, diagrams, 6 X 9in. Englewood 
Cliffs, New Jersey, Prentice-Hall, Inc., 
1960. Price: $12.00 (trade); $9.00 (text). 


This volume is the first of its kind to offer a 
complete block diagram analysis of each oper- 
ational unit of a digital computer—describing 
and illustrating in detail just how a computer 
is put together and how it works. The author 
emphasizes the organization and functional 
interrelation of fundamental units and dis- 
cusses these relationships from the viewpoint 
of operational necessity rather than mathe- 
matical logic, Boolean algebra or circuit 
theory. Especially helpful to the reader is 
the extensive treatment given two often- 
neglected topics: (1) input/output and (2) 
control of computer function. Problems are 
included that provide a means of acquiring 
facility with computer terms, and the final 
chapter reviews in detail what happens from 
the time the programmer receives the problem 
until the answer is ready. A complete glos- 
sary of terms is provided at the end of the 
book. 


ELECTRICAL ENGINEERING PROBLEMs, by E. 
M. Williams and F. J. Young. 294 pages, 
diagrams, 6 X 9 in. New York, McGraw- 
Hill Book Co., Inc., 1960. Price, $7.75. 


A collection of engineering problems ideally 
suited for laboratory activities is presented 
in this volume. Ranging from intermediate 
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to advanced the problems include such sub- 
jects as linear and nonlinear d-c and a-c cir- 
cuits, electromechanical devices, transformers, 
magnetic-core switching circuits, rotating- 
machines, and vacuum and semiconductor 
electronic circuits. The problems have all 
been laboratory tested during the past decade 
in various courses at Carnegie Institute of 
Technology which stress laboratory problem- 
solving disciplines. The book is particularly 
suitable for electrical engineering work, but 
can be used in the laboratory of a variety of 
courses, such as circuits, electronics, and 
mechanics. 


PROCEEDINGS OF THE SECOND CONFERENCE 
on REACTIONS BETWEEN COMPLEX NUCLEI, 
edited by A. Zucker, E. C. Halbert and 
F. T. Howard. 319 pages, illustrations, 
81 X11 in. New York, John Wiley & 
Sons, Inc., 1960. Price, $7.00. 


This book contains all forty research papers 
presented at this conference held in Gatlin- 
burg, Tennessee, May 2-4, 1960, and includes 
transcripts of the discussions which followed 
them. The papers are grouped into five 
sections to correspond to topics of the five 
conference sessions—transfer reactions; Cou- 
lomb excitation; scattering; transuranium, 
fission and low energy Li reactions; and com- 
pound nucleus reactions. In addition there 
is a summary of the conference by Prof. E. P. 
Wigner. These theoretical and experimental 
papers by distinguished scientists from Can- 
ada, Great Britain, Switzerland, USSR, as 
well as the United States, represent the results 
of various research projects, interpretation of 
these results, and discussion of their relation 
to the present body of knowledge about 
atomic nuclei. 
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Tuyratrons, by C. M. Swenne. 73 pages, 
plates, diagrams, 53 X 8} in. New York, 
The Macmillan Company, 1960. Price, 
$3.00. 


This book is written for the engineer who 
is already familiar with electrical engineering 
but who needs further study to keep up with 
industrial electronics. In four chapters the 
author covers physical principles, electrical 
characteristics, basic circuits, and applications 
of thyratrons (in relay circuits, electronic 
timers, controlled rectifiers, control circuits 
and inverters). The presentation is not 
cluttered up with mathematics; rather, it is 
a practical treatise, with descriptions and ap- 
plications set forth in plain language. An 
N. V. Phillips’ publication, it is published in 
English (Macmillan), French, German, Span- 
ish and Dutch. 


STATISTICAL THEORY AND METHODOLOGY IN 
SCIENCE AND ENGINEERING, by K. A. 
Brownlee. 570 pages, 6X9 in. New 
York, John Wiley & Sons, Inc., 1960. 
Price, $16.75. 


In preparing this elementary text on statis- 
tical methods, the author has “. . . leaned 
towards attempting to convey understanding 
of principles rather than practice, a choice 
rationalized on the ground that at this mathe- 
matical level the market is better provided 
with textbooks for the latter than the former.” 
The student is assumed to know exponents, 
logarithms, and simple algebraic manipulation. 
The main topics covered are discrete distribu- 
tion, nonparametric tests, analysis of variance, 
and regression and correlation. The Ap- 
pendix contains ten tables of useful data, in- 
cluding percentage points of the ¢, x2 and F 
distributions, random sampling numbers, etc.. 
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PUBLICATIONS RECEIVED 


FROM THEORY TO PRACTICE IN Sort MECHANICS. Selections from the writings of Karl Terzaghi, 
with bibliography and contributions on his life and achievements prepared by L. Bjerrum, 
A. Casagrande, R. B. Peck and A. W. Skempton. 425 pages, illustrations, 83 X 11 in. 
New York, John Wiley & Sons, Inc., 1960. Price, $12.00. 


Bounpary LAYER THEORY, by Hermann Schlichting, translated by J. Kesten. Fourth edition, 
647 pages, illustrations, 63 XK 9} in. New York, McGraw-Hill Book Co., Inc., 1960. 
Price, $16.50. 


HANDBOOK OF LAPLACE TRANSFORMATION: TABLES AND EXAMPLES, by Floyd E. Nixon. 
115 pages, 5 X 8 in. Englewood Cliffs, Prentice-Hall, Inc., 1960. Price: $6.00(trade) ; 
$4.50 (classroom). 


PERIODICALS IN THE CHEMICAL SocreTY LipraRyY. 48 pages, 6} X 8} in. London, The 
Chemical Society, 1960. Price, 5 shillings. 


Ciassics IN SocioLoGy. A course of selected reading by authorities, with an introductory 
reading guide by Donald MacRae. 326 pages, 53 X 9 in. New York, Philosophical 
Library, 1960. Price, $6.00. 


CLassics IN WESTERN CrviLizATIon. A course of selected reading by authorities, with an 
introductory reading guide by D. C. Somervell. 338 pages, 53 X 9 in. New York, 
Philosophical Library, 1960. Price, $6.00. 


CuLassics iN BrotoGy. A course of selected reading by authorities, with an introductory 
reading guide by S. Zuckerman. 351 pages, 5} X 9 in. New York, Philosophical 
Library, 1960. Price, $6.00. 


CLassics IN SCIENCE. A course of selected reading by authorities, with an introductory 
reading guide by E. N. Da C. Andrade. 322 pages, 5} X 9in. New York, Philosophical 
Library, 1960. Price, $6.00. 


READABLE RELATIVITY, by Clement V. Durell. 146 pages, diagrams, 5} X 8 in. New York, 


Harper & Brothers, 1960. Harper Torchbook No. 530. Price, $1.25 (paper). 


A History OF THE THEORIES OF AETHER & ELEctricITY, by Sir Edmund Whittaker. New 
York, Harper & Brothers, 1960. Harper Torchbooks Nos. 531 and 532. Vol. I: The 
Classical Theories, 434 pages, 5} X 8 in. Price, $1.95 (paper). Vol. Il: The Modern 
Theories, 319 pages, 53 X 8 in. Price, $1.85 (paper). 


CONCEPTS OF SPACE, by Max Jammer, with foreword by Albert Einstein. Reprinted from the 
1954 edition. 208 pages, 5} X 8 in. New York, Harper & Brothers, 1960. Harper 


Torchbook No. 533. Price, $1.40 (paper). 


THE FouNDATIONS OF ARITHMETIC, by Gottlob Frege, translated by J. L. Austin. 119 pages, 
53 X 8 in. New York, Harper & Brothers, 1960. Harper Torchbook No. 534. Price, 


$1.25 (paper). 


LECTURES ON THEORETICAL RHEOLOGY, by M. Reiner. Third edition, 158 pages, illustrations, 
53 X 8Lin. New York, Interscience Publishers, 1960. Price, $4.85. 
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THE FRANKLIN INSTITUTE 


Minutes of the Stated Meeting 
December 21, 1960 


The Stated Meeting of The Franklin Institute was held at 8:15 P.M. in the Lecture Hall, 
plewoaire by a dinner in Franklin Hall at which Silver Members of 1960 were honored. Wynn 

aurence LePage, President, called the meeting to order with approximately 120 members and 
guests present. 

The President stated that the minutes of the Stated Meeting of October 19, 1960 were 
published in the November issue of the JOURNAL. There being no corrections, the minutes 
were approved as published. He stated the minutes of the Stated Meeting of November 16, 
1960 will be published in the next issue of the JoURNAL and will be presented at the Annual 
Meeting in January 1961 for approval. 

The President announced that at the Annual Meetirg of The Institute to be held on 
January 18, 1961, ten members are to be elected to the Board of Managers for a term of three 
years beginning in 1961, and that the Secretary had informed him that eleven nominees had 
been named and nominations were closed fifteen days before the Stated Meeting in accordance 
with the By Laws. For the Chairman of the Nominating Committee, the Secretary read the 
following report: 

‘Under the provisions of Article VIII, Section 3 and Article IV, Section 5, of the By Laws, 
the following eleven names have been properly placed in nomination to fill the ten vacancies 
on the Board of Managers which will occur at the Annual Meeting of The Institute on January 
19, 1961. 


Morton Gibbons-Neff 1944* Henry B. Bryans 1951 Emery W. Loomis 1954 
Rupen Eksergian 1945 Wynn Laurence LePage 1953 Henry M. Chance, II 1958 
Francis J. Chesterman 1949 Gaylord P. Harnwell 1954 oseph A. Fisher 1959 
James Creese 1950 enneth Robert Atkins .... 


“The Nominating Committee is required by the By Laws to submit to the membership at 
this meeting its recommendations with respect to the candidates. 

“Accordingly, the Nominating Committee recommends to the membership that the follow- 
ing ten gentlemen be elected to fill the vacancies which will exist. All the recommended 
nominees are currently members of the Board of Managers: 


Morton Gibbons-Neff _— Creese Emery W. Loomis 
Rupen Eksergian enry B. Bryans Henry M. Chance, II 
Francis J. Chesterman Wynn Laurence LePage Joseph A. Fisher 


Gaylord P. Harnwell 


“The ballot will include all the nominees and the recommendations of this Committee will 
be indicated.” 

The President approinted the rong Tellers of Election, in accordance with the By Laws, 
who are neither officers nor nominees: Alfred W. Duerig, Dr. Joseph S. Hepburn, John H. Neher. 

Mr. LePage spoke briefly of some of the things that are in prospect for The Franklin 
Institute. Taking seriously the announcement of the President of the United States some 
eighteen months ago that science education was a national policy, and giving serious thought 
to the part The Institute has played in science education during the past 136 years, the Man- 
agers requested a survey be made to appraise the work of the Institute in this field. This was 
made possible by a grant from the 5 tmedi McMahon Foundation. Using this survey, a 
long-range program hee science education was developed. To give financial support to carry 
out this program, we have the assistance of leading citizens of the Philadelphia area, includin 
one of our own members, John A. Diemand, President of The Insurance Company of North 
America, who has consented to be General Chairman of the Drive for funds for The Franklin 
Institute Development Program. But The Franklin Institute is more than a science museum 
and planetarium. Here we conduct research and development work for industry and the 
government. We have launched a program to procure long-needed facilities for The Franklin 
Institute Laboratories for Research and Development and we will enter into a long-term lease 
for a research center to be built in an area adjacent to The Franklin Institute’s main building. 

In presenting the speaker of the evening, Miss Beth Peterson, Home Economist for the 
Extension Division of the E. I. duPont de Nemours & Company, Mr. LePage commented that 
he believed Miss Peterson was the first lady speaker to grace our platform. Miss Peterson’s 
talk on “Fibers, Facts and Fashions” was a fascinating description of the newer fibers and 
fabrics and the research required not only to manufacture the fabrics but to determine con- 
sumer requirements. Following the talk, members and guests were invited to inspect the 
samples of fabrics which were beautifully displayed in the front of the hall. 

The President expressed warm thanks to Miss Peterson for her delightfully interesting 
talk, and adjourned the meeting at 9:50 P.M. 

WILLIAM F. JAcKsON, JR. 
Secretary 


* Year first elected to Board of Managers. 
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New Books {J. F. I. 


NEW BOOKS IN THE FRANKLIN INSTITUTE LIBRARY 


MATHEMATICS 


FUNDAMENTALS OF SIGNAL THEORY, by J. L. Stewart. McGraw-Hill, 1960. 

OPERATIONS RESEARCH AND SYSTEMS ENGINEERING, by C. D. Flagel. Johns Hopkins Press, 
1960. 

SAMPLE DEsIGN IN BusINESS RESEARCH, by W. E. Deming. Wiley, 1960. 

THE THEORY OF STORAGE, by P. A. P. Moran. Wiley, 1959. 

ADVANCED ENGINEERING MATHEMATICS, by C. R. Wylie. 2nd ed. McGraw-Hill, 1960. 

E.Lectronic Computers, by T. E. Ivall. 2nded. Philosophical Library, 1960. 

OPERATIONS RESEARCH, by J. H. Batchelor. 2nd ed. Saint Louis University Press, 1959. 

INFORMATION AND DECISION Processes, by R. E. Machol. McGraw-Hill, 1960. 

DiGiTaL COMPUTERS AND NUCLEAR REACTOR CALCULATIONS, by W. C. Sangren. Wiley, 1960. 

CONTRIBUTIONS TO PROBABILITY AND STATISTICS; essays in honor of Harold Hotelling, edited 
by Ingram Olkin. Stanford University Press, 1960. 

INTRODUCTION TO LAPLACE TRANSFORMS FOR RADIO AND ELECTRONIC ENGINEERS, by W. 
D. Day. Interscience, 1960. 

STATIONARY PROCESSES AND PREDICTION THEORY, by H. Furstenberg. Princeton University 
Press, 1960. 

INTRODUCTION TO LINEAR PROGRAMMING, by W. W. Garvin. McGraw-Hill, 1960. 

RIEMANN SuRFACES, by L. V. Ahlfors. Princeton University Press, 1960. 

INTERNATIONAL CONFERENCE ON INFORMATION PROCESSING. Proceedings. UNESCO, 1959. 

DESIGN FOR A BRAIN, by W. R. Ashby. 2nded. Wiley, 1960. 

INTERNATIONAL CONFERENCE ON SCIENTIFIC INFORMATION. Proceedings. National Research 
Council, 1959. 

AUTOMATIC LANGUAGE TRANSLATION, by A. G. Oettinger. 


Harvard University Press, 1960. 


PHYSICS 


THE SCATTERING AND DIFFRACTION OF WAVES, by R. W. King. Harvard University Press, 
1959. 

INTRODUCTION TO QUANTUM MEcHANICs, by R. H. Dicke. Addison-Wesley, 1960. 

REACTOR ANALYsIS, by R. V. Meghreblian. McGraw-Hill, 1960. 

CONCEPTS OF THERMODYNAMICS, by E. F. Obert. McGraw-Hill, 1960. 

NucLear SPECTROSCOPY, edited by F. Ajzenberg-Selove. Academic Press, 1960. 

WAVE PROPAGATION IN A RANDOM Mept1uM, by L. A. Chernov. McGraw-Hill, 1960. 

ELEMENTS OF THE KINETIC THEORY OF GASES, by E. A. Guggenheim. Pergamon Press, 1960. 

BIBLIOGRAPHY ON PLASMA PHYSICS AND MAGNETOHYDRODYNAMICS AND THEIR APPLICATIONS 
TO CONTROLLED THERMONUCLEAR REACTIONS, compiled by J. D. Ramer. Maryland 
University, Engineering and Physical Sciences Library, 1959. 

Utrrasonics, by B. Carlin. 2nd ed. McGraw-Hill, 1960. 

CaREERS IN Puysics, by A. W. Smith. Rev. ed. Long’s College Book Co., 1960. 

Mecuanics, by K. R. Symon. 2nd ed. Addison-Wesley, 1960. 

Cryopnysics, by K. Mendelssohn. Interscience, 1960. 

BIOPHYSICAL SCIENCE, edited by J. L. Oncley. Wiley, 1959. 

PROPERTY MEASUREMENTS AT HIGH TEMPERATURES, by W. D. Kingery. Wiley, 1959. 

EINFUHRUNG IN DEN WARME UND STOFFAUSTAUSCH, by R. G. Eckert. 2. Aufl. Springer, 
1959. 

DIFFRACTION GRATINGS AS MEASURING SCALES, by J. Guild. Oxford University Press, 1960. 

THERMAL ENGINEERING, by H. L. Solberg. Wiley, 1960. 


CHEMISTRY AND CHEMICAL ENGINEERING 
ADVANCED ORGANIC CHEMISTRY, by G. W. Wheland. 3rded. Wiley, 1960. 
Liguips AND Liguip Mrxtures, by J. S. Rowlinson. Academic Press, 1959. 
PHYSICOCHEMICAL MEASUREMENTS AT HIGH TEMPERATURE, by J. O’M. Bockris. Academic 
Press, 1959. 
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THE SURFACE CHEMISTRY OF METALS AND SEMICONDUCTORS, edited by H. C. Gatos. Wiley, 
1960. 

EXPLORING THE STRUCTURE OF MATTER, by J. J. Trillat. Interscience, 1959. 

SYMPOSIUM ON THE MECHANISM OF HETEROGENEOUS CATALyYsIs. Proceedings. Elsevier, 
1960. 

PuLp AND Paper, by J. P. Casey. 2nded. Interscience, 1960. 

SUBSTITUTION AT ELEMENTS OTHER THAN CARBON, by C. K. Ingold. Weizmann Science 
Press of Israel, 1959. 

CHEMISTRY OF SYNTHETIC DruGs, by P. May. 5th ed. Longmans, 1959. 

POLYTHENE, by A. Renfrew. 2nd ed. Interscience, 1960. 

POLYSTYRENE, by W. C. Teach. Reinhold, 1960. 

MOLECULAR STRUCTURE, by J. C. D. Brand. Arnold, 1960. 


ELECTRICAL ENGINEERING 


TRAVELING-WAVE ENGINEERING, by R. K. Moore. McGraw-Hill, 1960. 

ELECTRONICS AND NUCLEONICS DICTIONARY, by N. M. Cooke. McGraw-Hill, 1960. 

A SIMPLIFIED TECHNIQUE OF CONTROL SYSTEM ENGINEERING, by G. K. Tucker. Minneapolis- 
Honeywell Regulator Co., Brown Instruments Division, 1960. 

PRINCIPLES OF FEEDBACK CONTROL, by C. H. Wilts. Addison-Wesley, 1960. 

ANALYSIS AND DESIGN OF FEEDBACK CONTROL SysTEMS, by G. J. Thaler. 2nded. McGraw- 
Hill, 1960. 

DicTIONARY OF AUTOMATIC CONTROL, by R. J. Bibbero. Reinhold, 1960. 


MATERIALS SCIENCE 


Cermets, by J. R. Tinklepaugh. Reinhold, 1960. 

Fiat RoLLep Propucts II, edited by E. W. Earhart. Interscience, 1960. 

Trin AnD Its ALLoys, by E. S. Hedges. Arnold, 1960. 

REPORT ON STRESS-CORROSION CRACKING OF AUSTENITIC CHROMIUM-NICKEL STAINLESS 
STEELS. American Society for Testing Materials, 1960. 

SYMPOSIUM ON PARTICLE SIZE MEASUREMENT. American Society for Testing Materials, 1959. 

SYMPOSIUM ON ELECTRON METALLOGRAPHY. American Society for Testing Materials, 1960. 

SYMPOSIUM ON EDUCATION IN MATERIALS. American Society for Testing Materials, 1960. 

SYMPOSIUM ON EFFECT OF WATER-REDUCING ADMIXTURES AND SET-RETARDING ADMIXTURES 
ON PROPERTIES OF CONCRETE. American Society for Testing Materials, 1960. 

SYMPOSIUM ON DURABILITY AND WEATHERING OF STRUCTURAL SANDWICH CONSTRUCTIONS. 
American Society for Testing Materials, 1960. 

SyMPosIuM ON NEWER METALs. American Society for Testing Materials, 1960. 

SYMPOSIUM ON FATIGUE OF AIRCRAFT STRUCTURES. American Society for Testing Materials, 
1960. 

SYMPOSIUM ON TREATED Woop FOR MARINE USE. American Society for Testing Materials, 
1960. 

DEsIGN oF STEEL STRUCTURES, by B. Bresler. Wiley, 1960. 

SYMPOSIUM ON THE THEORY OF THIN ELasTic SHELLS. Proceedings. Interscience, 1960. 

METALLIC CorROsION INuIBITORS, by I. N. Putilova. Pergamon Press, 1960. 

FLIESSPRESSEN VON STAHL, by H. D. Feldmann. Springer, 1959. 


4 
: 
2, 
3 
2 
we 


Artificial Reverberation.—|mprove- 
ments in artificial reverberation prom- 


ise to facilitate conversion of an 
ordinary auditorium into the acoustic 
equivalent of a full-sounding concert 
hall. This acoustic conversion, which 
introduces no undesirable side effects, 
is now being actively investigated at 
Bell Telephone Laboratories. As with 
other artificial reverberation, it is ac- 
complished entirely by electronic 
means. 

Called ‘‘colorless’’ artificial rever- 
beration, its principles were outlined 
at a recent meeting of the Audio 
Engineering Society. Dr. M. R. 
Schroeder and Mr. B. F. Logan of the 
Laboratories’ technical staff are jointly 


perfecting the new electroacoustic 
technique. 
Artificial reverberation has been 


used for many years by recording 
companies, broadcasters, and room 
acousticians. For example, a typical 
user might be a radio station wishing 
to broadcast a recording of an outdoor 
concert or open-air summer music 
festival. At the broadcast studio, an 
artificial reverberator would be em- 
ployed to introduce reverberations, or 
echoes, in an attempt to imitate the 
fullness of an indoor performance. 
The usual artificial reverberator 
available today has detrimental fluc- 
tuations in its amplitude-frequency- 
response characteristic which deviate 
appreciably from that of a real room. 
Thus an unpleasant quality may be 
added to the reverberant sound, an 
effect sometimes similar to talking 
through a long barrel. To obtain 


natural-sounding reverberation, users 
have had to resort to reverberation 
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chambers or other analog devices for 
producing the desired echoes. 

However, based on extensive studies 
of the problem, Schroeder and Logan 
of the Laboratories’ Visual and Acous- 
tic Research Department employed 
all-pass filters with exponentially 
decaying transients. The resulting ar- 
tificial reverberation responds equally 
to all frequencies. Spectral balance, 
or color of the reverberated sound 
is not altered. Hence it is called 
“colorless.” 

A concert-hall effect can be approxi- 
mated in the home through artificial 
reverberation. To achieve this with 
a radio or phonograph two or more 
loudspeakers are desirable. 

Messrs. Schroeder and Logan have 
conducted many colorless artificial re- 
verberation experiments at the Lab- 
oratories. To facilitate their research, 
they are investigating various schemes 
using the magnetic core storage and 
the logic units of a large-scale digital 
computer as the artificial reverberator. 
The reverberated sound produced as 
the computer’s output is evaluated 
subjectively by listening. 

A related subject also under study 
by Dr. Schroeder is quasi-stereophony. 
This method achieves the fullness (but 
not the localization) of true stereo- 
phonic sound with only one audio 
channel. Electroacoustic devices that 
produce this eftect without discoloring 
the sound have been successfully built 
and demonstrated at the Laboratories. 


A New Atomic Clock.—Three Har- 
vard scientists have built a new kind 
of ‘“‘maser’’ which promises to be over 
100,000 times more accurate than the 
best atomic clock now running. 
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The heart of the new device is an 
ingenious ‘bouncing box’’ that stores 
high-energy-level hydrogen atoms long 
enough for them to emit a detectable 
wave of radiation of precisely 21-cm. 
wavelength. 

The maser can be used both as a 
sensitive amplifier or as a source of 
21-cm. radiation of unprecedented sta- 
bility and purity of tone. The nearly 
perfect ‘‘monotone”’ it emits promises 
to be the most accurate timepiece in 
the world, a new standard for time 
or for, frequency. 

The accuracy of the new atomic 
maser has not yet been measured since 
with the best clock in the world, there 
is only one way to test its accuracy: 
build another identical clock and see 
how well they keep time with each 
other. A second atomic maser is now 
being built. It is hoped eventually 
that the accuracy will be one part 
in 10%. 

The new maser built by the Harvard 
scientists is the first to use free atoms, 
rather than atoms trapped in molecules 
or crystals. The radiation that the 
free hydrogen atoms give off is much 
purer than that emitted by other 
masers ; interaction among atoms held 
in a molecule or a crystal tends to blur 
or spread out the frequencies they emit. 

The hydrogen atom is made up of a 
proton and an electron. Because both 
are spinning, each acts as a tiny mag- 
net. In a hydrogen atom, the mag- 
netic axis of the proton can point 
either the same way as, or the opposite 
direction from, the magnetic axis of 
the electron. When the two magnetic 
axes flip from the parallel or high 
energy orientation to the anti-parallel 
position, the atom emits a photon of 
21-cm. wavelength. 

Here is how the atomic hydrogen 
maser works: A stream of hydrogen 
atoms is sent through a magnetic field 
that reflects those hydrogen atoms at 
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the high energy state into a unique 
storage box. Since the radiation given 
off by the hydrogen is very weak, the 
trick is to keep the atoms in the high 
energy state for a long time—until 
enough have accumulated to give off a 
detectable chorus of 21-cm. radiation. 

Ordinarily, the’ hydrogen atoms, 
bouncing off the sides of the box, 
would interact with the surface and 
quickly drop to the lower energy level. 
The new device has a special ‘‘bounc- 
ing box’’ a quartz bulb coated on the 
inside with a paraffin that does not 
very readily interact with the bounc- 
ing atoms. The scientists have found 
that one high-energy-level hydrogen 
atom can endure at least 10,000 bumps 
in this storage box, a life-time of about 
one second. 

In this way, a large number of high- 
energy-level hydrogen atoms is stored 
up. The gas will then spontaneously 
emit a nearly perfect ‘‘monotone,”’ 
thus acting as an oscillator or source 
of 21-cm. waves. But by deliberately 
decreasing the number of atoms pres- 
ent, this spontaneous emission may 
be checked. Then the device will 
emit 21-cm. radiation only when it is 
triggered by an incoming signal of the 
same wavelength, making the “‘bounc- 
ing box” a powerful amplifier. 

An atomic clock, essentially an 
atom or a molecule that emits radia- 
tion, is the most accurate kind of 
timepiece known. Its “‘pendulum” 
is the very steady stream of pulses, 
or waves, it sends out. 

As an oscillator, the new device 
will be used to measure some of the 
subtle physical properties of the pro- 
tons, electrons and neutrons in atoms 
much more accurately than has been 
possible before. The Harvard men 
plan to study the interaction between 
the nucleus and its orbiting electron 
in deuterium and tritium (heavy iso- 
topes of hydrogen). 
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In another experiment, they will 
deliberately force into the storage box 
an excess of hydrogen atoms to study 
and measure the interactions among 
them. These same __ interactions 
occur in the great clouds of hydrogen, 
in and between galaxies, that also 
broadcast on the 21-cm. wavelength. 

The device could also be used to 
measure, to a much higher precision 
than has yet been possible, the effect 
of both gravity and speed on time, 
which was predicted by Einstein. 

The atomic maser might eventually 
be of use in long-range communication. 
Such an amplifier in an earth satellite 
could be used to bounce signals from 
earth to another point. Because the 
amplifier’s radiation is so pure, its 
range would be enormous. 

The research was sponsored by the 
National Science Foundation and the 
Office of Naval Research. 


Plasma Accelerator Tested.—A 
“plasma gun,” built to test new ways 
of keeping future space vehicles travel- 
ling steady during flight, was fired 
continuously nearly four million times 
during a test at General Electric’s 
Missile and Space Vehicle Department 
in Philadelphia. 

Scientists at the department’s Space 
Sciences Laboratory said the machine 
operated without stop for 184 hr. and 
estimated it provided enough thrust 
during the test to stabilize a far-trav- 
elling spacecraft for two years. The 
thrust, or recoil, from the firings would 
be used to stop unwanted spinning 
and tumbling motions. 

The equipment is called REPPAC 
which stands for Repetitively Pulsed 
Plasma Accelerator, and its test for 
more than 18 hours is believed by the 
laboratory to have set a record for 
this type of apparatus. 

The new type of plasma propulsion 
will have a big advantage over the 
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more conventional rockets and gas 
jets which have already been success- 
fully used to stabilize space vehicles 
and satellites, according to General 
Electric scientists because : 

The plasma accelerator shoots off a 
very energetic stream of material at 
extremely high speeds to obtain thrust, 
and a small amount of this material 
(plasma) goes a long way. Conven- 
tional gas jets or rockets obtain thrust 
by shooting off comparatively large 
amounts of material at much lower 
velocities, and thus much more fuel 
or compressed gas has to be stored. 
As the length of time of the space 
flight increases the jet and rocket 
methods become more undesirable 
because of weight and storage space 
limitations. 

The plasma (a gas whose atoms 
have been broken up) reacts with a 
powerful magnetic field which acceler- 
ates the gas particles 50 times a second 
and discharges them at more than 
22,000 miles per hour to produce 
thrust. 

The REPPAC’s “gun” is a T- 
shaped tube in which small amounts of 
gas are turned into the plasma and 
propelled into a vacuum chamber by 
a combination of electrical charges 
and a magnetic field. 


All-Glass Insulating Paper.—New 
applications are demonstrating the 
ability of AMF’s Tissuglas—an all- 
glass insulating paper—to solve un- 
usual electrical, thermal, and electronic 
problems. Tissuglas, a product of 
American Machine & Foundry Com- 
pany, is available in thicknesses as 
low as 0.0006 of an inch. It is uniform 
in texture and electrical resistance. 

The newest, and potentially largest, 
application of Tissuglas is as a base 
for electroluminescent lamps. These 
lamps are flexible panels as thin as 7 
of an inch. The resulting lamps are 


= 
j 
i 
x 
3 : 
q 
ae 
: 
a 
i 


Jan., 1961.] 


said to be the brightest, clearest and 
most uniformly textured light of any 
electroluminescent panels on the mar- 
ket today. 

Another new application of Tissu- 
glas is in the manufacture of bonded 
strain gages. Tissuglas, in this appli- 
cation, is used as a saturating base. 
In this use, its ability to pick up 
resins in uniform thicknesses (a func- 
tion of its patented floc-free matting) 
gives the gages uniform resistance and 
high temperature stability. 

Tissuglas is also finding a new mar- 
ket as a component for printed circuit 
materials. In this application, Tissu- 
glas is laminated with a very thin cop- 
per foil. As an example, the one mil 
Tissuglas together with the ? mil cop- 
per produce the thinnest and most 
flexible printed circuit material now 
available. 

Tissuglas is a sheet of matted sub- 
micron glass fibers made on a specially 
adapted paper making machine. It 
has very unusual properties as an in- 
sulating material: melting tempera- 
ture above 1,200 degrees F.; thick- 
nesses from ultrathin (.0006 in.) to 
.012 in.; pore sizes from 8 to 90 
microns; and a dielectric of 370 to 
2,370 volts. Produced in any width 
up to 38 in. and in continuous rolls. 


Electronic Telephone Service.—A 
customer trial of the world’s first elec- 
tronic telephone central office (ECO) 
was started in mid-November of 1960 
by Bell Telephone Laboratories and 
Illinois Bell Telephone, at Morris, 
Ill. 

Developed by Bell Laboratories, it 
is the first system which performs 
electronically all the functions in 
handling of telephone calls. It was 
described by Bell Laboratories Presi- 
dent James B. Fisk as“ . . . a totally 
new system making possible a variety 
of new services which will double the 
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flexibility and hence the potential 
value of each customer’s telephone.” 

Through the use of tiny gas tubes, 
a vast memory system, and thousands 
of transistors and other solid state 
devices, which operate in millionths 
of a second, ECO permits Morris trial 
users to: 

—Use home extension phones as 
intercoms, simply by dialing two 
digits. 

—Reach frequently called numbers 
by dialing only two digits instead of 
seven. 

—Have incoming calls routed to 
another phone when the original called 
line is busy. 

—Dial a code which causes all sub- 
sequent incoming calls to be auto- 
matically transferred to any other 
number. (For example, to a friend’s 
home if you went there for the 
evening.) 

Other new services expected to be 
introduced after the trial is under way 
will permit many of the trial cus- 
tomers to: 

—Allow a third party to be called in- 
to an existing telephone conversation. 

—Have an immediate connection 
to a busy line as soon as the line be- 
comes available. 

Information gained from the Morris 
trial will tell not only how ECO per- 
forms, but it will also give an insight 
into which of the new service features 
telephone users prefer. Telephone 
engineers are now investigating ways 
to provide some of these new services 
through existing telephone systems. 

The trial will also present an oppor- 
tunity to study the kinds of mainte- 
nance activity and the number of people 
needed to operate and maintain the 
new exchanges. During the trial, the 
Morris ECO will be operated jointly 
by Illinois Bell and Bell Laboratories. 

In appearance, ECO has little in 
common with present exchanges. The 
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long, high, racks of equipment and 
the noisy chatter of relays are gone. 
Instead, there are several rows of neat 
gray cabinets filled with thousands of 
plug-in packages of electronic com- 
ponents. The switching network, 
made up of 23,000 tiny, neon-filled 
tubes, goes about its work quietly, 
with only a bright red glow to show 
which tubes are operating. 

One of the keys to the success of 
the ECO is that it uses stored program 
control to perform its job—that is, 
each action of the machine is deter- 
mined by instructions stored in its 
memory. 

The difference between ECO’s 
stored program control and present 
systems can be simply illustrated by 
the two different ways a student can 
solve a square root problem. One way 
is by a series of logical steps—pointing 
off two places, etc. The other way is 
by looking up the answer in a book of 
square root tables. 

Existing exchanges use the first 
method to work out the detailed logic 
of each telephone call, item by item. 
The ECO, on the other hand, solves 
many of its problems by “looking up 
the answer in the book.” In this case 
the “‘book”’ is information—2} million 
bits—stored on photographic plates 
in the machine’s memory, telling it 
what to do in any situation. 

Because of ECO’s extreme flexi- 
bility, it is capable of offering new and 
different services. 

The ECO performs its work with 
such speed that it has enough “‘spare”’ 
time to continually check its own 
circuits. When it discovers a fault, 
it locates and diagnoses the trouble 
and, in some cases, even makes the 
necessary corrections. If the fault 
can’t be corrected by the machine, it 
writes out in detail on a Teletype- 
writer what the trouble is, adding the 
month, day, hour, and minute of the 
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malfunction. Then the repairs are 
made by a member of ECO staff. 

Although some of its basic theory 
has been known for many years, ECO 
became feasible only with the inven- 
tion of the transistor and other semi- 
conductor devices, and with the de- 
velopment of high speed, high capacity 
memory systems. To illustrate the 
importance of semiconductors to the 
system it has only to be noted that 
the Morris ECO uses some 12,000 
transistors and 105,000 diodes. It 
uses other electronic components by 
the tens of thousands. 

Although the Morris ECO is an 
experimental version, Bell engineers 
are already at work on a production 
model, suitable for volume manufac- 
ture by the Western Electric Com- 
pany. Many components in the pro- 
duction model will differ from those 
used at Morris, but the over-all philos- 
ophy of the two systems will be simi- 
lar. The first production model will 
be operating by mid 1965, and from 
that time on ECO will gradually be 
introduced into the Bell System to 
meet the needs of telephone growth 
and to replace worn-out central offices. 


Teaching Devices Demonstrate 
Physics Principles.—Two new labora- 
tory demonstration pieces suggested 
by a prominent instructor to help 
teach physics students about changes 
created by thermal energy are now 
being made by Central Scientific Co., 
Chicago. One of the devices demon- 
strates the forces involved in thermal 
expansion and contraction while the 
other shows the effects of temperature 
change on the resistance of conductors. 

Both were developed by Julius Sum- 
ner Miller, professor of physics at El 
Camino College, California, who 
brought them to Central Scientific. 
Cenco makes and markets many simi- 
lar laboratory instruction pieces. 
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The thermal expansion and contrac- 
tion unit is the No. 77430 Cenco- 
Miller Forces in Thermal Expansion 
Apparatus. This unit is built around 
a steel expansion bar about 1 ft. long, 
held between two ends of a yoke. 
The rod can be heated by a manifold 
gas burner built into the apparatus; 
it can be cooled either by the sur- 
rounding air or by spraying water 
on it. 

In the test, cast iron drive pins are 
inserted in the rod at the free end, 
one on each side of the yoke end piece. 
As the heated rod expands, its pres- 
sure is exerted against the inside drive 
pin, causing it to break with an audible 
snap. On contraction, the other pin 
will snap. 

Built of heavy duty material for 
years of use, this large demonstration 
(163 in. long) can be seen throughout 
the classroom. 

The other demonstration is the No. 
83064 Cenco-Miller Temperature Co- 
efficient of Resistance Apparatus. In 
this device, one of the two wires lead- 
ing to a 6-volt 32-candlepower lamp 
consists of resistance wire wound 
around two rods so that the wire forms 
a bed. A tubular gas burner beneath 
this wire bed heats it, causing the 
resistance to increase and the light 
to go out. Cooling causes the lamp 
to brighten. This unit is 16} in. long. 

In addition to these demonstrations 
of thermal effects, two other Cenco- 
Miller devices demonstrating rota- 
tional forces are now available from 
Cenco. For additional information 
about any of the new Cenco-Miller 
demonstration devices, write to Cen- 
tral Scientific Company, 1700 Irving 
Park Road, Chicago 13, Illinois. 


Scrambling Device Blocks Bank 
Passbook Forgery.—A scrambling de- 
vice which makes it virtually impos- 
sible to forge passbook signatures in 
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the withdrawal of savings bank de- 
posits and, at the same time, sub- 
stantially reduces customer waiting 
time at tellers’ windows has been an- 
nounced by the Radio Corporation of 
America. 

Known as Signaguard, the device 
reproduces a passbook signature as 
an unrecognizable mass of broken 
lines. When the passbook is presented 
at the teller’s window, the device re- 
turns the signature to its original ap- 
pearance for comparison with the cus- 
tomer’s signature on the withdrawal 
slip. 

Dr. H. J. Wall, Manager of RCA 
Applied Research at Camden, N. J., 
which developed Signaguard, said 
banks are expected to use the device 
in conjunction with electronic data- 
processing systems situated in their 
main offices. 

This would eliminate the need for 
duplicate records of signature and ac- 
count status in each branch office, he 
pointed out, and also would do away 
with time-consuming withdrawal 
checking procedures. 

The new protection system, Dr. 
Wall said, makes use of fiber optics— 
glass tubes that carry light and images 
around bends and corners. 

It so effectively scrambles a signa- 
ture that even a bonafide depositor, 
much less an aspiring forger, cannot 
decipher the name on a lost passbook, 
he said. The imprint is produced as 
mixed segments of lines thousandths 
of an inch in diameter, corresponding 
to the diameter of the tubes. 

As a double check, a secret number 
can be attached to the name. This 
effectively foils forgery attempts by 
persons who know the depositor well 
and can duplicate his signature. 

The new device, Dr. Wall said, will 
make savings bank accounting fool- 
proof. Another of its advantages is 
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taining cards on all depositors—which 
each teller must consult before com- 
pleting a withdrawal. This, he said, 
would make for more room and less 
confusion in the banks, and would 
eliminate the time required to refer 
to such files for signature comparison 
and account checking. 

Signaguard may be likened to a 
cable or telephone wires, with each 
individual wire being a glass tube in- 
stead of copper wire, the RCA official 
explained. Each tube picks up a small 
segment of the signature and trans- 
mits it to the other end of the tube 
via a devious route. This scatters sig- 
nature segments throughout the un- 
intelligible mass that is imprinted on 
sensitized paper at the other end of 
the tube. 

In the bank the fiber optics tube is 
reversed to bring the scrambled signa- 
ture back to its original form for com- 
parison with the signature submitted 
to the teller. 

Use of Signaguard in conjunction 
with the bank’s central computer sys- 
tem, Mr. Wall said, would permit a 
teller to complete a banking transac- 
tion without leaving his window. 

Applications, other than in savings 
banks, are contemplated for Signa- 
guard, Dr. Wall said. It may be em- 
ployed in identification cards for com- 
mercial bank depositors, for retail 
credit cards, and in-plant security. 


Thermocouples for Low Tempera- 
ture Applications.—A new line of 
gold-cobalt vs. copper thermocouples 
and complementary indicating, record- 
ing and controlling instruments for 
low temperature applications has been 
announced by Cryogenics, Inc., Wash- 
ington, D. C. 

Open and closed well thermocouples 
are available in ranges from 0-300°K, 
with accuracy of +1 per cent of span 
or 2°K, whichever is greater. All are 
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interchangeable to +4°K. In special 
cases, the company notes, a spool of 
one continuous length prepared from 
the same melt is provided and held for 
customers, with all thermocouples be- 
ing fabricated from this spool. In 
such situations, interchangeability of 
better than 2°K is guaranteed. 

Direct Reading Instruments with 
internal reference junction compensa- 
tion for use with these thermocouples 
are also available, in a thoice of indi- 
cating, recording or _ controlling 
functions. 

For additional information, write 
Cyrogenics, Inc., 1129 Vermont Ave., 
N. W., Washington 5, D. C. 


Commercial Gold Recovery Unit.— 
The first commercial unit for recover- 
ing gold lost in electroplating has been 
introduced by Technic, Inc., of Provi- 
dence, R. I1., one of the largest pro- 
ducers of precious metal solutions for 
plating. 

Until now between 4 and 10 per 
cent of all gold used in electroplating 
has simply gone down the drain, and 
the loss to the plating industry runs 
to millions of dollars a year. 

With the new unit, one electroplater, 
using the new special resin Gold Saver 
on a trial basis, will have recovered 
enough gold in less than three months 
to pay for the unit. 

The Gold Saver, designed to become 
part of the gold-plating process, con- 
sists basically of : 

—A column of special resin, small 
“BB” size particles, which capture 
gold in solution that now is being lost ; 

—A pump for circulating through 
the resin column the water in which 
the gold is dissolved ; and 

—Pipes, hoses, and connections be- 
tween the pump and the special resin 
column. 

After an item is plated with gold, 
it is immediately immersed in a tank 
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of water to rinse off and save the 
plating solution, which frequently 
costs between $40 and $50 a gallon. 
This “rinse” tank is called a drag-out 
tank. 

Over a period of time, the concen- 
tration of gold solution in the drag-out 
tank builds up until, because of its 
gold content, it is worth possibly $10 
a gallon. Periodically contents of this 
tank are added to the gold plating 
tank to keep the plating solution at a 
workable level and to recover some 
of the gold. 

In the past, the plated item was 
then rinsed in running water, and the 
gold “drag-out” from the drag-out 
tank was lost. 

The Gold Saver comes before this 
running water rinse. 

The unit is attached to a tank into 
which the plated item is dipped, after 
the drag-out tank. The pump con- 
tinuously circulates the solution 
through the special resin trap, where 
the gold is picked up. 

The work, free of gold, is then rinsed 
in running water. 

When the resin is completely satur- 
ated with gold it may be burned by 
the plater to recover the gold, sent to 
a refiner, who would also recover the 
gold by burning, or it may be returned 
to the Technic laboratories for 
reclaiming. 

The Gold Saver is designed for any 
acid or non-free cyanide golds, such 
as Technic’s 24K acid, Orotherm HT, 
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Orosene 999, or HG non-free cyanide 
gold baths. Similar ion-exchange 
traps for cyanide golds, rhodium, pal- 
ladium, and platinum already are 
under development or test. 


New Cancer Treatment.—A new 
method for using cytotoxic drugs in 
the treatment of cancer has been de- 
veloped by surgeons at the Royal In- 
firmary, Glasgow, Scotland. It has 
been applied successfully in four 
operations. 

The use of cytotoxic drugs, de- 
veloped in the U. S. as a cancer pal- 
liative by reducing the size of tumors, 
or in conjunction with surgery, has 
been limited up to now because it can 
be lethal: a dose strong enough to kill 
a tumor can sometimes also kill the 
patient. 

The Glasgow surgeons have de- 
veloped a technique to eliminate this 
danger by isolating the circulation in 
the affected part of the body. By so 
doing, 20 to 30 times the amount of 
cytotoxin can be administered than 
previously. 

The surgeons are achieving pallia- 
tion for two to four years of cancers 
which were previously incurable, and 
they are also able for the first time to 
investigate the properties of the cyto- 
toxic substance. A hospital official 
said, ‘“‘The work is not being dupli- 
cated elsewhere.” 

The method is still in the develop- 
ment stage. 
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YOU CAN ADVANCE SCIENCE EDUCATION 


Today, more than ever before in its 135-year history, there is vital need for 
The Franklin Institute to effectively promote education in science and technology. 
It is imperative that we meet this challenge by providing adequate educational op- 
portunities in these fields. This requires vision, objective planning, and money. 
We have more than enough of the first two requisites, but far too little of the third. 

Our programs are aimed at professional scientists and industry, as well as the lay 
public and young people seeking inspiration and guidance in choosing a career. 
The Institute’s educational programs are impressive, for they begin with students in 
the early grades of our elementary schools and continue throughout an individual’s 
professional or industrial life. With more funds at our disposal, the scope and 
vigor of these activities could be greatly increased. 

The Franklin Institute is not richly endowed. It is a non-profit organization, 
depending for encouragement and support on an understanding public. Capable 
and conservative management assures wise administration of all funds. 

Your gift or bequest, large or small, will be deeply appreciated and will be used 
effectively to broaden the Institute’s educational usefulness. There is a warm 
satisfaction in giving financial support to an organization that has pioneered in, and 
is dedicated to, the advancement of science and technology through education. 

When property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

The Secretary of The Franklin Institute will gladly furnish you with additional 
information. Write to him at The Franklin Institute, Benjamin Franklin Parkway 
at Twentieth Street, Philadelphia 3, Pennsylvania. 
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THE FRANKLIN INSTITUTE 
LABORATORIES FOR RESEARCH AND DEVELOPMENT 


announce opportunities for: 


ELECTRICAL ENGINEERS 
SYSTEMS ENGINEERS 
ORGANIC CHEMISTS 
OPERATIONS RESEARCH SPECIALISTS 
MECHANICAL ENGINEERS 
EXPERIMENTAL PSYCHOLOGISTS 
METALLURGISTS 
INDUSTRIAL ENGINEERS 
PHYSICISTS 
APPLIED MATHEMATICIANS 


to work on challenging problems in research 


Send complete resume to: 
Mr. John E. Christ, Director of Personnel 


THE FRANKLIN INSTITUTE 
Philadelphia 3, Pa. 
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BUSINESS 


MATH ALY 
ENGINEERING 


Now in the fourth year of operation, The Franklin Institute Com- 
puting Center has shared its technical know-how, modern electronic 
equipment and experienced personnel with hundreds of progressive 
industries and government agencies across the nation. 


We offer the services of creative people, highly skilled in their respec- 
tive fields and ably trained in the application of these skills to the ever 
expanding area of electronic computers and data processing systems. 


This staff is now available for analysis, system design, programming 
or coding of projects of unlimited scope or context. Input to our large 
scale computer and peripheral equipment is acceptable in any form. 
Results are provided on cards, plastic or metallic tape, and in com- 
pletely edited printed copy. Our extensive library of automatic coding, 
engineering, data processing and mathematical routines is available to 
all users, and machine time is provided with or without the services 
of programming personnel. 


THE FRANKLIN INSTITUTE 
Computing Center 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pa. LOcust 4-3600 
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